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I.  INTRODUCTION 


The  purpose  of  this  thesis  is  to  determine  an  optimum 
Naval  Airborne  Tactical  Electronic  Warfare  System  from 
selected  components.  Its  characteristics  of  performance  will 
be  derived  from  particular  mission  requirements,  generated 
from  Ecuadorian  naval  tactical  operational  needs. 

In  history  there  are  plenty  of  examples  that  confirm 
that  a  military  action  will  better  succeed  if  sufficient 
reconnaissance/surveillance  has  supported  the  final  decision. 
In  this  particular  case,  an  ESM  system  must  provide  sufficient 
and  accurate  information  in  order  to  locate  the  enemy  and 
clarify  the  tactical  situation,  that  is;  what  are  the  enemy 
forces  disposition,  strength,  intentions,  discretion, 
readiness,  etc.?.  In  other  words,  the  ESM  receiver  should 
collect  data  from  all  the  emitters  of  interest,  over  a  wide 
frequency  spectrum,  over  large  angular  regions,  and  with 
sufficient  sensitivity  to  improve  detection  ranges;  then  it 
must  measure  and  process  this  information  quickly  and 
accurately  for  immediate  use  as  a  decision  basis. 

There  are  basically  three  different  kinds  of  receivers 
available  in  a  very  mature  technological  state:  a.)  Crystal 
Video  Receivers  (CVR),  b.)  Instantaneous  Frequency 
Measurement  (IFM)  and  c.  )  Superheterodyne.  The  particular 
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characteristics  of  each  one  of  them,  will  be  analyzed  to 
identify  their  applicability  to  Naval  requirements. 

A.  BACKGROUND 

In  the  past  twenty  years,  there  has  been  a  growing 
interest  among  Third-World  nations  in  strengthening  the  naval 
presence  in  oceanic  areas  of  interest,  justified  by  principles 
and  policies  of  military  strategy,  and  due  to  the  natural 
growth  or  maturity  achieved  with  time  and  experience. 

The  late  sixties  brought  about  challenging  new  weapons: 
the  Sea  skimming  surface  to  surface  missile  (SSM’s)  ,  and  the 
fast  attack  missile  boat  (FAMB’s).  As  a  result  the  concepts 
of  Naval  Warfare  were  shaken  with  the  enormous  changes  that 
these  two  weapons  represented  to  the  traditional  way  of  making 
war  at  sea. 

Subsequently,  both  the  FAMB’s  and  the  SSM’s  were  proven 
in  different  actions  around  the  world,  and  due  to  their 
successes  many  Navies  have  adopted  them  as  the  "pattern  to 
follow",  even  though  in  many  cases  they  were  not  wisely 
utilized,  or  were  not  used  for  those  purposes  for  which  both 
the  missile  or  the  FAMB  were  appropriate.  There  are  many 
examples  to  cite,  where  countries  have  acquired  FAMB’s  for 
purposes  of  patrol  ling  their  Exclusive  Economic  Zone  but  found 
them  to  be  not  cost-effective.  In  other  cases,  endurance 


and  seakeeping 

capabi 1 i ti es 

of 

the 

FAMBs 

have  been 

overestimated . 

The  fact  is 

that 

the 

presence 

of  FAMB’s 
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represents  a  threat  to  forces  attempting  to  penetrate  their 
defensive  zones.  Additionally,  their  cost  of  operation  is 
much  less  than  for  a  Destroyer  or  Frigate. 

Ecuador  is  a  particular  case  where  FAMB’s  were  bought  in 
order  to  replace  Second  World  War  vintage  destroyers, 
achieving  at  low  cost  a  considerable  power  projection 
capability.  This  was  not  the  best  solution  for  two  reasons: 

First  because  of  the  geographical  situation  of  that 
country,  (located  at  the  Northwestern  part  of  South  America 
and  with  a  wide  open  sea  in  front  of  it)  was  not  the  best 
environment  for  missile  boat  operations  since  they  are  better 
used  in  restricted  waters.  And  second,  the  large  open  ocean 
area,  combined  with  the  limited  endurance  and  surveillance 
capabilities  of  the  FAMB’s,  made  it  very  impractical  to 
search,  locate,  attack  and  destroy  threats  with  the  required 
level  of  success.  As  with  bigger  ships,  the  problem  of 
searching  and  locating  the  enemy  remains  the  key  point  for 
success  at  sea,  and  this  is  very  difficult  to  soive  with  the 
limitations  of  the  FAMB’s.  In  any  case,  during  the  UNITAS 
operations  held  from  1984  to  1986  (in  which  the  author 
participated),  an  experimental  solution  was  operationally 
tested  which  consisted  of  combining  a  Surface  Attack  Group 
with  a  P-3C,  where  the  aircraft  carried  the  task  of  searching 
for  and  locating  enemy  forces  and  then  vectoring  the  missile 
boats  towards  the  threat  until  they  reached  missile  launching 
position. 
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It  may  be  important  to  note  that  operating  with  this 
combination  gave  very  significant  results  in  extending  the 
effective  operational  radius  of  the  FAMB’s  and  their 
lethality,  especially  when  engaging  surface  forces  with 
limited  air  support.  Clearly,  the  aircraft  tasks  were  of 
extreme  importance,  and  particularly  the  ESM  tactics  and 
techniques  employed  which  were  paramount  to  overall  mission 
success.  In  addition  to  this,  the  aircraft  enjoyed  the 
advantage  of  its  versatility  which  allowed  it  to  selectively 
operate  its  active  emitters,  and  thereby  approach  closer  to 
the  target  (under  certain  conditions)  for  positive 
identification. 

It  is  becoming  evident  that  a  cost-effective  solution  for 
improving  ones  own  power  projection  capabilities  is  to  outfit 
surveillance  aircraft  with  appropriate  ESM  systems,  and  to 
develop  the  related  tactics  for  their  employment  either  alone 
or  in  combination  with  a  Surface  Attack  Group.  This  thesis 
will  analyze  several  airborne  ESM  System  in  order  to 
satisfactorily  accomplish  the  mission  and  operational 
requirements  as  stated  above. 
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I  1 .  SCENAR 1 0 


The  scenario  considered  by  this  thesis  is  that  of  a  Naval 
aircraft  in  Electronic  Warfare  surveillance  operations 
attempting  to  target  an  already  classified  enemy  ship.  The 
aircraft  will  have  the  support  of  an  Surface  Attack  Group 
(SAG)  for  weapons  delivery  on  the  target. 

The  object  of  the  ESM  mission  is  to  obtain  and  process 
sufficient  and  accurate  ESM  data.  This  is  not  a  simple  matter 
since  many  complex  parameters  are  involved.  Some  parameters 
may  no  be  improved  without  trading  off  the  performance  of  the 
system  in  another  area.  Other  parameters  are  under  the 
enemy’s  control  or  are  of  random  nature.  This  thesis  is 
limited  to  a  particular  scenario,  which  is  divided  into  three 
specific  parts: 

-  The  Geographical  region, 

-  Ones  own  forces  and, 

-  The  possible  threat 

A.  THE  GEOGRAPHICAL  REGION 

The  Geographical  region  of  interest  is  the  Northwestern 
part  of  South  America,  (where  Ecuador  is  located)  from  the 
coastline  up  to  about  900  miles  westward,  and  Northward  to 
the  Panama  Canal.  Historically,  the  extent  of  the  region  has 
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been  justified  by  both  national  and  hemispherical  interests. 
Actually,  Ecuador’s  national  maritime  interests  are  the 
protection  of  its  Economic  Exclusive  Zone  and  maritime 
communication  lines;  and  in  the  case  of  an  hemispherical 


threat,  the  navies 

of 

this  region 

wi  1  1 

have  to 

share  the 

responsibi 1 ities 

for 

protection 

of 

vital 

maritime 

communication  lines.  In  fact,  the  most  recent  UN1TAS 
operations  have  underscored  this  point  by  increasing  the  joint 
exercises  between  the  navies  of  Ecuador,  Colombia,  Venezuela, 
and  the  USA,  on  both  the  Atlantic  and  the  Pacific  sides  of 
the  Panama  Canal.  Ecuador  possesses  naval  and  aeronaval 
facilities  along  the  Pacific  coast  and  in  the  Galapagos 
Islands  which  are  located  600  miles  west  of  the  continent. 

B.  FRIENDLY  FORCES 

The  available  forces  will  be  mostly  Missile  Boats 
ranging  upwards  from  125  Tons  to  800  Tons,  equipped  with  ESM 
equipment  on  board.  Actually,  Ecuador  possesses  a  missile 
force  of  6  Corvettes  and  6  fast  attack  missile  boats.  Since 
each  different  sized  FAMB  will  have  different  endurance,  it 
is  necessary  to  deploy  them  wisely  according  to  mission 
characteristics,  scale  of  conflict  and  nature  of  the  threat. 
In  any  case,  during  the  crucial  initial  phases  of  the 
encounter,  both  surveillance  and  command  and  control 
requirements  must  be  accomplished  by  the  ESM  aircraft  in  order 
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to  ach ; eve  the  required  degrees  of  surprise  and  accuracy  in 
the  surface  attack. 

The  main  weapon  chosen  was  the  French  Exocet  MM38  and 
MM40  missiles  with  21  and  40  nautical  miles  of  range 
respective ly. 

As  was  stated  above,  the  deployment  of  forces  depends 
on  complex  own  and  enemy  related  factors.  The  following  are 
examples  of  scenarios  that  might  happen  and  how  own  forces 
would  be  arranged,  according  to  range  from  the  advanced  bases 
avai lable. 

1.  Close  to  Coast  Engagements  (up  to  100  na. ) 

It  would  be  advisable  to  use  mostly  smaller  FAMB’s, 
with  the  size  of  the  threat  determining  how  many  boats  will 
be  dispatched  to  achieve  a  power  advantage. 

2.  Medium  Range  Engagements  (lOO  to  200  na. ) 

In  this  scenario,  the  general  choice  is  to  select 
larger  FAMBs.  However  if  enemy  classification  is  obtained  in 
time,  it  might  be  advisable  to  send  a  combination  of  FAMB’s, 
with  the  smaller  units  closer  to  the  target,  thereby  allowing 
simultaneously  firing  of  the  MM38  and  MM40  missiles.  With 
this  deployment,  one  can  exploit  the  proven  advantages  of 
relative  Radar  and  ESM  detection  ranges,  using  the  smaller 
FAMB’s  with  their  reduced  radar  cross  sections  to  proceed 
closer  to  their  targets. 
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3.  Long  Range  Engagements  (>  200  nm.  ) 

Longer  endurance  units  (Corvettes)  will  be  required 
for  this  scenario,  but  the  attack  group  must  also  include  a 
sufficient  number  of  FAMB’s  to  counter  the  threat  size. 

Some  scenario  parameters  evaluated  as  common  to  all 
the  cases  described  above  ares 

a)  The  surveillance  aircraft  must  detect,  classify, 
identify  and  locate  the  enemy  ship  prior  to  further  action; 
and  additionally  to  this,  it  must  vector  and  command  the 
attack  group  properly  until  the  point  where  Boats’  commanders 
assume  total  control  of  the  action. 

b)  A  minimum  of  2  attacking  units  will  be  required. 

c)  A  minimum  of  2  missiles  will  be  fired  at  each 
target  to  increase  the  probability  of  hit. 

d)  After  the  attack  has  been  made,  the  unit's  will 
leave  the  area  at  full  speed, 

e)  The  surveillance  aircraft  will  be  used  for 
evaluating  the  results  of  the  attack,  and/or  collecting 
valuable  information  that  might  support  conclusions  about  the 
outcome  of  the  action  or  previous  situation.  Using  this 
information  it  will  allow: 

1)  Reengagement  if  the  situation  is  favorable. 

2)  Relocation  of  additional  ships  if  saturation 
of  enemy’s  defenses  is  required, 

oi  . 
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3)  Launch  of  a  second  attack  if  the  enemy’s 
reaction  has  been  weak  and  slow, 
f)  Prior  to  and  during  the  attack,  the  deployed 
SAG  will  maintain  a  distance  between  units  of  about  5  to  15 
miles  for  ESM  triangulation,  UHF  communication  purposes,  and 
for  fire  power  coverage. 

All  these  courses  of  action  are  derived  from  the 
information  obtained  by  electronic  surveillance  during  the 
search  and  detection  phase  of  the  engagement. 

C.  THE  POSSIBLE  THREAT 

The  threat  which  this  thesis  analyzes,  is  a  surface 
warship  such  of  a  destroyer  or  cruiser  type,  acting  as  a 
picket  of  a  surface  force,  or  as  part  of  a  surface  action 
group,  and  having  little  or  no  airborne  early  warning  support. 
Additionally,  the  ship  will  have  an  on-board  he  1 o  for  anti- 
surface  operations. 

The  main  weap.  ns  system  for  this  ship  consists  of  a 
medium  range  AA  missile  (30  to  40  nm.),  a  short  range  AA 
missile  (10  nm.  )  and  over-the-horizon  surface  to  surface 
mi ss  i  1 es . 

In  order  to  properly  accomplish  these  tasks,  this  ship 
will  have  the  following  sensors: 

1  Long  range  air  search  radar  (frequency  agility 
capab 1 e ) , 
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1  Medium  range-low  altitude  air  search  radar, 
1  Tactical  surface  search  radar, 

1  Navigation  radar, 

3  Fire  control  radars,  and 
1  Surface  search  radar  on  the  helo. 
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III.  ESM  AND  THE  PARAMETERS  TO  BE  MEASURED 


A.  ENVIRONMENTAL  ELECTRONIC  DENSITY 

The  environmental  electronic  density  plays  a  key  role  in 
the  design  of  the  ESM  system,  since  this  density  establishes 
some  of  the  operational  requirements  under  which  the  system 
will  be  deve 1  oped . 

For  example,  the  number  of  emitters  and  their 
transmitting  characteristics,  their  location  relative  to  own 
receiver,  the  level  of  hostilities,  and  the  type  of  mission 
to  be  accomplished  by  own  forces,  are  among  all  the  factors 
that  outline  the  environment  one  will  work  with. 

According  to  Ref.C5]  the  electromagnetic  density  for 
systems  with  sensitivity  below  -70  dBm  is  particularly 
altitude  independent.  This  particular  value  deserves  to  be 
mentioned  because  we  will  be  refer  back  to  it  later  when  the 
ESM  parameters  are  defined  and  analyzed  within  each  System 
conf i guration. 

From  the  model  cited  above,  and  using  the  formula  for 
maximum  unambiguous  range,  we  can  estimate  the  pulse  density 
to  be  the  following  (pulses  per  second): 

PRF  =  c/2  •  range  (2-1) 
where : 
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PRF 


=  Pulse  repetition  Frequency  (pulses  per  second) 
c  =  Speed  of  1 i ght 
range=  Maximum  Unambiguous  Range 


The  radars  expected  to  be  present  are: 


a ) 

Long  range  search  radar  (aerial) 

300 

Fp. 

(277nm. range) 

b) 

Medium  range  search  (aerial) 

1012 

Fp 

(80 

nm. range ) 

c ) 

Tactical  surface  search 

1079 

Fp. 

(75 

nm. range ) 

d) 

Navigation  radar 

1620 

Fp. 

(50 

nm  range) 

e ) 

Fire  control  radars  (3) 

2699 

Fp. 

(30 

nm  range) 

f  ) 

Search  radar  (Helo) 

1000 

Fp. 

(80 

nm  range; 

This  makes  a  total  of  13,108  Fp  for  each  ship  in  the  case 
where  all  the  sensors  are  radiating.  This  number  give  us  a 
good  starting  point  for  the  analysis  since,  as  stated  above, 
we  might  expect  to  find  more  than  one  enemy  ship  in  the  zone, 
in  addition,  merchant  ships  will  be  in  the  area  and  at  some 
stage  we  will  also  have  to  count  own  ships  radars.  Assuming 
that  we  will  encounter  a  maximum  of  10  ships  with  similar 
characteristics  to  those  cited  above,  all  transmitting  their 
radars  simultaneously,  we  will  have  a  total  Pulses  per  Second 
of  131,080  Fp  (pps). 

For  the  purposes  of  this  thesis  the  airborne  system  will 
have  a  processing  capability  of  about  300,000  pulses  per 
second,  in  order  io  satisfy  and  exceed,  within  a  sufficient 
level,  the  requirements  imposed  by  the  tactical  situation. 


12 


B.  PARAMETERS  OF  INTEREST 

In  order  to  create  a  decision  basis,  an  ESM  system  must 
detect  and  process  many  parameters  from  the  emitter,  and  the 
more  parameters  are  obtained  and  successfully  analyzed,  the 
easier  is  to  determine  the  emitter's  identity. 

Additionally,  some  ESM  parameters  are  considered  more 
valuable  than  others  in  a  given  situation;  such  as  the  Angle 
of  Arrival  (AOA)  or  the  frequency  characteristics  of  the 
intercepted  signal,  which  contribute  the  most  to  emitter 
identification.  Generally,  the  ESM  receiver  will  be  required 
to  measure  the  following  parameters  from  the  emitter: 

Angle  of  arrival  (AOA), 

Frequency  ( f  )  , 

Time  of  arrival  (TOA), 

Pulse  width  (PW),  and 
Pulse  amplitude  (PA). 

These  parameters  have  been  intentionally  ordered 
according  to  their  importance  to  mission  accomplishment 
(i.e.,  sorting  and  identification). 

As  can  be  found  on  the  commercial  market,  ESM  receivers 
range  from  detecting  one  to  detecting  all  of  the  above 
parameters,  but  if  emitter  sorting,  classification  and 
identification  are  required  (which  is  usually  the  case  with 
naval  ESM)  a  minimum  of  three  parameters  are  needed. 

A  basic  receiver  scheme  is  independent  of  which 
parameters  are  being  measures.  For  example: 
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Fig.  1,  Basic  receiver  system 


C,  DYNAMIC  RANGE  AND  SENSITIVITY 

Dynamic  Range  and  Sensitivity  are  of  extreme  importance 
in  ESM  receiver  design.  If  performance  in  both  are  poorly 
realized,  the  overall  system  performance  will  be  seriously 


degraded. 

On 

the  other  hand, 

if  good 

dynamic 

range 

and 

sens i t i vi ty 

are 

avai 1 ab 1 e  at 

the  receiver, 

then 

many 

parameters 

such 

as  detection 

range, 

accuracy  of 

the 

measurements,  sorting  capability,  probability  of  detection, 
etc.  will  give  the  system  and  the  user  an  advantageous 
decision  basis.  In  the  author’s  point  of  view,  any  design 
must  first  start  by  setting  appropriate  values  for  these  two 
fundamental  specifications. 


14 


1.  Dynamic  Range 

The  Dynamic  Range  defines  the  range  of  input  power 
over  which  the  receiver  will  work  properly,  generally 
expressed  in  decibels  (dB).  The  lower  limit  of  the  dynamic 
range  is  the  system’s  sensitivity,  and  the  upper  is  defined 
in  different  manners,  the  most  common  being  the  one  dB 
compression  point.  The  input  power  is  linearly  related  to  the 
output  power,  but  as  the  input  power  increases  there  is  a 
certain  limit  where  the  output  power  loses  its  linear 
relationship  and  falls  below  the  ideal.  When  the  difference 
between  the  ideal  linear  receiver  and  the  actual  output  power 
is  1  dB,  then  the  value  of  the  input  power  (at  which  the  1  dB 
difference  occurs)  is  the  upper  limit  of  the  required  dynamic 
range. 

If  the  input  power  goes  above  this  limit,  the 
receiver  saturates,  the  gain  is  lowered,  and  spurious  signals 
will  appear  at  the  Receiver’s  output.  To  cope  with  the 
expected  range  of  the  signals  strengths,  a  Dynamic  Range  of 
at  least  65  dB  is  required. 

2.  Receiver  Sensitivity 

The  System’s  Sensitivity  defines  the  weakest  signal 
power  that  can  be  detected  by  the  System  and  generate  an 
output  where  the  Signal  to  Noise  ratio  meets  the  minimum 
requ i red . 
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As  will  be  addressed  later  in  the  detectability 
problem,  a  certain  minimum  power  level  of  the  target  emitter 
is  also  required  to  produce  a  detection  at  a  given 
probability.  The  emitter’s  radiation  power  must  compete 
inside  the  ESM  receiver  with  the  intrinsic  thermal  noise 
(internally  generated)  that  represents  the  floor  of  the 
system.  The  sensitivity  floor  is  a  function  of  both  noise 
power  and  signal  to  noise  ratio  (SNR). 

The  noise  power  (N)  at  the  receiver’s  output  is 
modeled  by  the  well  known  following  relationship  that  assumes 
a  Gaussian  noise  (resistor)  model: 


N  =  k • To ■ Bn • Fn 


(3-1  ) 


where : 


-23 

k  =  1.38  •  10  Boltzman’s  constant 

To  =  Noise  Temperature  (in  Kelvin) 

Bn  =  Receiver  noise  Bandwidth  (Hz.) 

Fn  =  Receiver  noise  Figure 


If  we  use  T  =  290°  K,  and  express  the  power  in 
mill i watts  then : 


N(dBm)  =  -174  dBm  +  10  log,*  (Bn)  +  10  log,.,  (Fn)(3-2) 

10  10 


The  receiver  noise  figure  (Fn)  relates  the  actual 
noise  power  output  to  the  ideal  receiver’s  output . 
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The  effective  noise  bandwidth  (Bn)  becomes  a  ver.y 
important  consideration  when  instantaneous  wideband  receivers 
are  used,  therefore,  even  though  a  wide  bandwidth  might  be 
available  for  interception  at  the  receiver’s  front  end,  this 
does  not  represent  directly  the  noise  power  available  because 
when  the  signal  reaches  the  detector,  most  of  the  original 
noise  is  filtered  out  by  the  narrow  bandwidth  of  the  video 
stage.  This  is  referred  to  sometimes  as  "Processing  Gain", 
and  it  is  proportional  to  the  ratio  of  the  Acceptance 
Bandwidth  (Ba)  to  the  Video  Bandwidth  (Bv),  where  gamma  is  the 
integration  efficiency. 


Bn 


Brf 


Brf 


2  Bv 


gamma 


(3-3) 


Gamma  varies  from  0.5,  (when  the  ratio  Brf/Bv  is 
large)  to  about  0.8  or  0.9  for  Superheterodyne  receivers 
(  Brf/Bv  =  2,  approx.). 

The  signal  received  from  the  emitter  is  given  by: 
Pt*Gt*Gesm*\^»g^ 

S  =  -  (3-4) 

2 

(4*tc»R)  •  Lesm*Lt 


where: 


Pt  =  Radar  transmitter  peak  power  (watts) 

Gt  =  Radar  transmitter  antenna  gain 

Gesm  =  ESM  receiver  antenna  gain 

\  =  c/f;  Wavelength  (meters) 

g  =  Multipath  factor 

R  =  Range  between  radar  and  ESM  receiver 
Lesm  =  ESM  receiver  losses 
Lt  =  Radar  transmitter  losses 

The  Signal  to  noise  ratio  will  be  defined  by 
combining  the  previous  three  equations: 

Pt*Gt*Gesm»\2*g2 

SNR  =  -  (3-5) 

2 

(4kR)  •k»To»Bn»Fn*Lt»Lesm 

Where  Bn  is  in  Hertz. 

These  two  definitions,  noise  floor  and  SNR,  are 
quite  important  because  they  represent  the  basis  for  the 
following  system  design.  Two  parameters  directly  related  to 
the  noise  and  SNR  of  the  system  are  the  maximum  detection 
range  and  the  probability  of  detection  /  false  alarm. 

The  power  advantage  of  the  ESM  receiver  comes  from 
its  1/R2  propagation  loss  versus  the  1/R^  propagation  loss 
for  the  radar. 
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In  Fig.  1  we  see  the  relationship  between  signal 
power  received  at  the  ESM  and  radar  receivers  in  dBw  as  a 
function  of  range,  showing  a  difference  that  increases  with 
distance. 

In  appendix  A,  MATHCAD  2.0  is  used  for  solving  the 
equations  and  plotting  the  variations  of  signal  power  at 
reception,  both  at  the  Radar  and  at  the  ESM  receiver,  given 
their  characteristics.  Figure  2  shows  the  1/1?  advantage  of 
the  ESM  receiver  which  allows  it  (given  the  same  range  for 
Radar  and  ESM  Receiver)  the  use  of  a  lower  sensitivity  level; 
and  if  the  sensitivity  of  the  ESM  receiver  is  enhanced  from 
this  value,  a  tremendous  advantage  in  detection  range  is 
achieved.  The  Radar  and  ESM  parameters  involved  in  the 
following  plots  are: 


Radar;  Long  Range 
Pt  =  250  Kw. 

Gt  =  1585 

PU  =  4.0  10"6  s. 

\  =  0. 1  m 

Lt  =  2 

g  =  1 

ht  =  20  m 


Air  Surveillance  type, 

Transmitter  Peak  Power  (Kw. ) 
Antenna  Gain  (mainlobe) 
Pulse  Width 
Wavelength  (meters) 

Transmitter  losses 
Multipath  factor 
Antenna  height 


ESM  System: 

Gesm  =  1  ESM  receiver  Antenna  Gain  (Omni) 
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hesm 


1000  m 


ESM  aircraft  height 


POWER 


Fig.  2  Radar  and  ESM  receiver  propagation,  in  dB,  for  a 
particular  radar. 

The  following  step  is  to  add  more  realistic  data  to 
the  comparison  by  introducing  the  radar  range  equation  to  the 
effects  of  noise  and  the  minimum  signal  to  noise  ratio 
required.  For  a  specified  radar  design,  the  most  important 


improvement  is 

achieved  if 

the  SNR  is 

set 

to  the 

1  owest 

possible  level, 

therefore, 

if  the  SNR 

i  s 

def i ned 

as  the 

independent  variable  of  the  equation,  we  can  estimate  the 
radar  detection  range.  In  this  analysis,  the  SNR  (minimum) 
of  the  radar  is  varied  from  10  dB  to  15.5  dB  (10  to  35.5) 
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radar  detection  range.  In  this  analysis,  the  SNR  (minimum) 
of  the  radar  is  varied  from  10  dB  to  15.5  dB  (10  to  35.5) 
which  are  values  commonly  employed  by  radar  designers. 

As  indicated  in  Figure  3,  using  a  conservative 
approach,  the  target  has  a  10  dB  SNR  at  a  maximum  detection 
range  of  about  200  Km..  This  suggests  why  it  is  important  to 
gather  basic  data  from  a  hostile  emitter  (peak  power,  antenna 
size  and  gain,  pulse  width,  frequency,  etc.  ),  since  it  will 
allow  E.W.  planners  to  estimate  radar  performance  and  counter 
with  the  appropriate  ESIi  equipment. 


Fig.  3  Radar  maximum  detection  range  (R  )  as  function  of 
°  rmax 

minimum  signal  to  noise  ratio  (SNR  .  ) 

6  min 
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detection  range  is  set  as  a  function  of  the  sensitivity  and 
the  sensitivity  is  varied  over  the  range  stated,  the  plot 
shown  in  Figure  4  is  obtained. 

From  Figure  4  it  can  be  concluded  that  for  Sensiti¬ 
vities  above  -45  dBm,  the  ESM  receiver  performance  is  poor 
since  we  know  that  the  maximum  detection  range  of  the  Radar 
is  about  200  Km.;  but  as  the  sensitivity  drops  to  -50  dBm  or 
below,  the  increment  in  ESM  detection  range  increases  rapidly. 
The  best  weapon  that  airborne  naval  ESM  surveillance  could 
make  use  of  is  its  detection  range  advantage.  Therefore,  no 
receiver  with  sensitivities  above  -50  dBm  will  be  considered 
appropriate  for  the  purposes  of  this  thesis  and  that 
automatically  eliminates  the  use  of  CVRs  for  further 
evaluaion. 


Fig  4.  ESM  performance  in  Range  for  detection 


D.  GEOMETRICAL  CONSIDERATIONS 

In  developing  a  tactical  ESM  mission,  it  is  important  to 
ensure  that  flight  altitude  and  detection  range  are  defined 
together.  To  determine  this  relationship,  there  are  three 
parameters  to  be  considered: 

a)  First,  given  the  emitter  and  aircraft  altitudes,  a 
radar  horizon  exists  that  physically  limits  microwave 
propagation.  The  following  equation  calculates  the  horizon 
limited  range,  dependent  on  Earth’s  sphericity  and  atmospheric 
refraction.  This  relationship  is  taken  from  [Ref.  43: 


R,  =  130.34  •  <  hesm  >’5  +  (  ht  )'5 

hoz  | 

where : 

ht  =  Radar  height  (in  Km) 

hesm  =  ESM  receiver  height  (in  Km) 


(3-6) 


b)  Second,  we  define  a  maximum  ESM  detection  range 

(R  )  dependent  on  ESM  receiver  and  radar  characteristics: 
esm 


0.  5 

Pt-Gt-Gesm-\2-g2  ] 


(4*n)2*S  .  -Lt-Lesm 
min 


(3-7) 


If  the  aircraft’s  altitude  is  too  low,  the  radar  horizon 
becomes  the  limit  for  the  maximum  detection  range. 
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in  order  to  optimize  the  capability  for  searching  and 
detecting,  the  flight  height  must  be  chosen  such  that  the 
radar  horizon  equals  the  maximum  ESM  detection  ranges 


(3-8) 


130. 34. 


0.  5 


(h 


esm 


(Km))'5+  (ht(Km))5 


2  2 

Pt- .  Gt-  *Gesm  •  \  •  g 


(4-rt)  *5  .  -Lt-Lesm 
min 


(3-9) 

h  ( Km ) 

esm 


r  _r  Pt*Gt*Gesm*  \  »g  1 

V  .  O 

7.67.10-6r 

-  ht(Km)0,5 

l  L(4*7t)*S  .  *Lt*Lesm- 

min 

■ 

This  is  a  basis  for  later  reference,  and  is  simplified, 
not  taking  into  account  other  factors  that  affect  the 
performance  of  the  system.  This  approach  assumes  the  maximum 
detection  range  as  the  principal  factor  of  interest  but  most 
of  the  time  it  is  limited  by  the  aircraft  altitude  (about 
10  Km.  ) . 

As  can  be  seen  in  the  Figures  5  and  6  which  plot  the 
height  for  the  radar  side  and  main  lobe  detection,  the 
altitudes  required  are  relatively  low  for  sensitivities  above 
-65  dBm  (300  m.  for  -65  dBm)  against  the  side  lobes,  on  the 
other  hand  if  the  main  lobe  gain  is  considered  in  the 
calculation,  the  height  goes  far  above  reasonable  or 
operationally  realistic  numbers  (above  1000  Km.  for  -75  dBm, 


and  around  13  Km  for  -65  dBm). 


c)  Final ly  a  most  important  consideration  is  grazing 
angle.  If  an  aircraft  is  located  at  or  below  the  tangential 
line  of  sight  defined  by  the  radar  horizon,  the  performance 
of  the  radar  emitter  will  be  considerably  degraded  because  of 
the  low  grazing  angle  of  the  emission  over  the  earth's 
surface.  This  increases  the  sea  clutter  and  diffraction 
effects  that  can  produce  a  lot  of  scattering  and  fading  on  the 
radar  signal  during  its  propagation  and  reception.  The  use 
of  this  tangential  line  will  define  what  will  be  called  the 
ESM  horizon  slant  range. 


h 

f  Igt 


3. 16 -10  min 

(-105  dBw ;  -75  dBm) 


1.010  watts 
(-80  dBw;  -50  dBm) 


Fig.  5  Flight  height  for  Main  lobe  detection 


Additionally  because  of 

the 

round 

tr  i"p 

requirement  it 

becomes  very  difficult  for 

the 

Radar 

to 

perform  in  an 

acceptable  way  under  such  conditions.  On  the  other  hand,  the 
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ESM  receiver  is  not  seriously  affected  by  this  phenomena 
because  it  will  receive  the  direct  ray  without  problems  and 


«  - - 

Km. 

10 

h 

[s  1 

f  lgt2 

L  min  J 
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V 

i, 

l 
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-12 

S 

-10 

10 

min 

3. 16  ■  10 

watts 

(-120 

dB; 

-90  dBm) 

(-95  dB.;  -65 

dBm) 

Fig.  6  Flight  height  for  side  lobes  detection 


with  a  relatively  high  power  level.  Ideally  it  would  be 
preferred  to  always  fly  in  this  tangential  approach.  Then, 
once  detection  has  occurred  while  flying  towards  the  emitter, 
it  must  continuously  lower  its  altitude  trying  to  stay  within 
the  radar  worst  detection  zone  the  maximum  possible  time. 
Within  the  scope  of  this  thesis,  the  aircraft  will  be  assumed 
to  be  flying  at  constant  speed  and  altitude,  since  the  purpose 
is  to  evaluate  and  analyze  the  variations  of  time  and 
reception  performance  parameters  in  a  simulated  environment. 
However,  further  research  should  include  the  tangential  flight 
path  combined  with  of  other  ESM  requirements  such  as  signal 


to  noise  ratio,  or  time  of  coincidences,  probability  to 
intercept  vs.  time,  etc.,  to  investigate  additional 
performance  improvements. 

In  Figures  7  and  8  we  see  the  theoretical  conditions  of 
this  tangential  approach,  which  includes  heights,  ranges,  main 
lobe  and  side  lobes  detection. 


Fig.  7  Electronic  Horizon  slant  range. 

We  can  conclude  at  this  stage,  that  the  performance  of 
an  airborne  ESM  system  initially  depends  on  how  far  away  it 
is  able  to  search,  therefore  the  flight  height  must  be 
properly  chosen  in  order  to  exploit  the  tangential 
approximation  the  longest  time  possible.  Lowering  the  range 
improves  the  emitter  side  lobes  interception  by  the  ESM 


receiver  but  the  altitude  must  be  adjusted  so  that  the  smaller 


range,  indicated  from  the  sensitivity  standpoint,  is  just  at 
the  electronic  horizon  slant  range.  If  the  altitude  is  higher 
at  the  time  this  side  lobes  detection  is  achieved,  we  might 
well  be  inside  the  radar  detection  range  in  a  region  of 
better  radar  reception,  having  lost  the  detectability 
advantage  of  the  ESM/Elint  aircraft. 

A  sensitive  analysis,  (Appendix  J)  using  MATHCAD 
demonstrated  that  the  improvement  of  radar  sidelobes  detection 
is  not  as  important  as  the  ability  to  obtain  the  maximum 
information  possible  from  the  main  lobe  emission  and  this  is 
a  strong  function  of  the  maximum  detection  range. 


POWER 


Fig.  8  Power  availability  at  radar  and  ESM  Rx.  versus  range 
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For  the  model,  depicted  in  Figure  7,  which  represents 
a  typical  radar  and  an  IFM  receiver,  detection  of  the 
sidelobes  (which  are  32  dB  below  the  main  lobe),  occurs  at 
ranges  less  than  the  radar  maximum  detection  range.  Perhaps 
a  more  sensitive  ESM  receiver  should  be  considered  if  sensing 
the  side  lobes  while  avoiding  radar  detections  is  such  a 
desirable  objective.  An  improvement  is  obtainable  with  their 
higher  sensitivity  receivers,  and  it  may  be  justifiable  for 
an  ESM  system  designer  to  evaluate  this  capability 
quant i tat i ve 1 y . 


E.  PARAMETERS  ACCURACY 

The  real  accuracy  of  the  ESM  parameters  of  interest  is 
obtained  by  statistical  analysis  of  sample  measurements, 
gathered  experimentally,  and  is  expressed  in  standard 
deviation  form. 


E  (  x.  -  x  ) 

l 


2  , 


0.  5 


(3-10) 


where : 


n  =  Number  of  measurements  taken 

x  -  parameter  Sample  mean 
x  . =  Sample  data  value 
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To  improve  the  accuracy  of  the  measurements,  n  should  be 
large  (  n>30  ).  The  quality  of  the  ESM  receiver  system  can 
be  inferred  if  the  difference  term  of  the  numerator  is  small 
and  if  the  parameter  sample  mean  is  close  to  the  source 
parameters . 

According  to  Wiley  [Ref.  2,  p.135]  the  accuracy  value, 
required  for  an  acceptable  measurement,  ranges  from  0.  IX  to 
IX,  which  is  the  ratio  of  the  standard  deviation  to  the  actual 
value  of  the  parameter  being  measured.  Ortiz  [Ref.  9],  and 
Giaquinto  [Ref. 10]  in  their  respective  works  used  accuracy 
requirement  values  of  0.  IX  for  frequency  measurement  and  IX 
for  both  Pulse  Repetition  Interval  (PRI)  and  Pulse  Width  (PW) 
measurements.  These  were  reasonably  chosen  and  correspond 
to  values  in  commercially  available  receivers. 

From  Refs.  3,  A,  9,  10  and  12,  the  theoretical  standard 
deviation  errors  for  frequency,  PRI,  and  PW  (as  functions  of 
the  Signal  to  Noise  Ratio)  for  a  Superhet  receiver  are 
respect i ve 1 y : 


df  = 


where 


3 

_  2*n2-PW3-Bn  SNR 
Ba  is  the  receiver 


0.  5 

;  max 


acceptance 


Ba 


bandwidth. 


(3-11;  3-12) 


dPR  I 


0  5 

0.35- (2) 

0. 8-Bv- (SNR)0, 5 


(3-13) 
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dPW 


0.7 


(3-14) 


Bv  [  2 • SNR  ]  °* 5 


The  formulas  above  described  are  derived  from  the 
analysis  done  by  Skolnik  [Ref.  4]  in  chapter  11  of  his  book; 
but  after  comparing  the  obtainable  standard  deviation  errors 
by  using  the  formulas  for  frequency  resolution,  It  was  found 
that  they  gave  very  different  values  which  need  to  be 
explained  for  application  to  an  ESM  receiver.  These  formulas 
apply  to  matched  receiver  superhets.  Normally  the  measurement 
accuracy  of  a  scanning  (non-mathed)  superhet  is  its  bandwidth. 
For  example,  a  superhet  tunable  over  a  range  (D)  of  one  GHz. 
and  having  an  acceptance  bandwidth  (Ba)  of  one  MHz.  would  have 
a  frequency  resolution  of  0.1%  of  its  tuning  range.  If  the 
tuning  range  were,  for  example  two  to  four  GHz.;  then  one 
could  have  a  0.1%  accuracy  with  a  bandwidth  of  two  MHz. 

For  IFM  receivers,  especially  the  digital  types, 
frequency  resolution  depends  on  the  number  of  discriminators 
or  the  least  significant  bit  (LSB)  of  the  multi-discriminator, 
binary  weighted,  digital  IFM. 

Operational  values  for  frequency  accuracy  in  IFM 
receivers  achieved  by  the  author  have  been  around  6  to  7  MHz 
and  that  is  also  mentioned  as  a  standard  frequency  accuracy 
by  C.  L.  Davies  in  his  article  [Ref.  14,  p.  165-1663  where 
he  talks  about  frequency  measurements.  In  general  terms  then, 
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the  real  accuracy  should  be,  as  stated  initially,  around  0.1 
to  1%  of  the  frequency  value. 

In  wideband  ESM  receivers  of  the  non-digital  IFM  type 
such  as  CVR  or  wide  band  superhets  -such  as  microscan 
receivers-  the  noise  bandwidth  will  be  considerably  less  than 
the  RF  bandwidth.  A  good  approximation  for  Bn  is: 

0.5 

Bn  =  [  2  -  B  •  Bv  ]  (3-16) 

In  any  case,  frequency  resolution  is  a  matter  of 
bandwidth  in  tunable  receivers.  In  IFMs  it  varies  from  very 
poor  (equal  to  Ba  in  CVR)  to  good  (0.1%)  for  digital  IFMs. 

The  last  parameter  to  be  considered  is  the  Angle  of 
Arrival  deviation,  (dAOA).  This  is  function  of  the  kind  of 
direction  finding  receiver  used  by  the  system,  and  for  the 
purposes  of  this  study  it  is  assumed  that  this  is  an 
independent  Receiver.  Actually  due  to  tne  wide  ESM 
requirements  in  angle,  power  and  frequency  (and  due  to  the 
platform  physical  limitations),  an  AOA  accuracy  of  5  to  7 
degrees  is  considered  to  be  typical  and  operationally 
reasonable  [Ref.  143. 
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IV.  ESM  MEASURES  OF  EFFECTIVENESS 


The  tactical  objectives  imposed  on  an  ESM  System  are  to 
search  for,  intercept,  locate  and  identify  all  the  microwave 
emissions  expected  to  be  found  in  the  probable  theater  of 
operations.  If  these  are  the  requirements,  a  Measure  of 
Effectiveness  (MOE)  must  reflect  quantitatively  (and  to  a 
certain  degree  of  accuracy),  how  well  an  ESM  System  is 
achieving  these  objectives. 

The  Probability  of  Intercept  (POI)  as  defined  by  Ortiz 
and  Giaquinto  [Refs.  9  and  103  is  a  function  of  three 
independent  factors:  detectability,  coincidence  and 
identification.  If  reasonably  good  values  are  obtained  for 
these  three  factors,  the  system  should  be  accomplishing  the 
tactical  requirements,  therefore  P01  appears  to  be  a  good  ESM 
system  measure  of  effectiveness.  The  model  for  the  POI  is  in 
fact  more  realistic  and  complete  than  most  of  the  definitions 
found  in  manufacturers  brochures  or  EW  texts.  It  emphasizes 
that  all  three  conditions  must  be  considered  if  a  sensitive 
analysis  is  needed  and  it  is  intrinsically  suitable  as  a 
decision  tool.  However,  it  is  important  to  keep  in  mind  that 
the  interception  problem  will  not  be  solved  by  just  having  a 
wide  reception  bandwidth,  high  sensitivity,  and  an  omni 
antenna:  it  is  necessary  to  quantify  the  relationships  and 
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trade  offs  between  them  and  then  identify  the  improvable  areas 
within  the  ESM  receiver. 

The  nature  of  the  variables  involved  in  this  problem  are 
classified  as  controllable,  uncontrollable  (but  predictable) 
and  random.  They  will  be  identified  as  the  model  is 
developed  with  the  fundamental  purpose  of  highlighting  their 
influence  on  the  final  output  and  how  an  ESM  designer  can 
overcome  the  induced  handicaps  by  using  his  controllable 
parameters . 


A.  PROBABILITY  OF  DETECTION 

The  first  independent  factor  considered  in  the  PO I  model 
is  the  probability  of  detection.  It  is  modeled  by  Skolnik 
tRef.  43  and  expanded  for  EW  applications  by  Tsui  [Ref.  3) 
where  the  ratio  of  IF  bandwidth  to  video  bandwidth  produces 
a  significant  improvement,  since  the  detector  behaves  as  a 
filter.  The  probability  of  detection  is  a  function  of  the 
signal  power,  receiver  sensitivity,  and  the  threshold  level, 
although  it  is  important  to  note  that  the  ESM  designer  can 
obtain  improvements  in  all  three  factors  by  devoting  attention 
to  the  variables  that  are  under  his  control. 

The  threshold  level  is  obtained  from: 


Pfa  =  exp 

where : 


=  Noise  voltage  variance 
Pfa  =  Probability  of  False  Alarm 


(4-1 ) 
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Generally,  the  probability  of  false  alarm  (Pfa)  is  set 
to  a  reasonable  value,  and  in  this  thesis  it  is  set  to  10 
Additionally,  the  noise  variance  is  set  to  1  [Ref.  3,  p.223. 

The  probability  of  detection  set  of  formulas  according 
to  Tsui  [Ref.  3,  Ch  2]  is  found  in  Appendices  D  and  E  of  this 
thesis  for  Superhet  and  1FM  receivers  respectively. 

It  should  be  noted  that  Tsui  has  introduced  the  ratio  of 
Br  to  Bv  in  the  probability  of  detection  equations,  therefore 
taking  in  to  account  what  is  called  as  the  "processing  gain" 
of  the  receiver.  Many  curves  are  shown  in  Tsui  [Ref.  3, 

Ch  23  where  the  overall  probability  of  detection  improvement, 
due  to  the  ratio  Br/Bv  (  T  in  the  calculations),  is 
emphasized.  The  available  signal  power  at  the  receiver,  and 
the  receiver  sensitivity,  utilized  in  Tsui’s  equations  were 
analyzed  in  Chapter  3  of  this  thesis. 

B.  PROBABILITY  OF  COINCIDENCE 

The  second  independent  factor  of  the  POI  is  the 
probability  of  coincidence. 

Before  a  detection  can  take  place,  the  receiver  must 
first  obtain  an  alignment  of  different  parametric  windows. 
Self  [Ref.  153  developed  a  simple  but  effective  model  that 
represents  each  parameter  involved  in  the  interception  process 
as  pulse  train  functions  of  time.  Each  window  function  is 
periodic,  will  have  a  window  width  (it  will  be  of  "squared" 
shape),  and  the  phase  or  starting  time  is  assumed  to  be 
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uniformly  distributed.  By  using  this  model  we  can  represent 
any  of  the  periodic  window  function  involved  in  the  problem 
such  as: 


-  Antenna  rotation  and  beamwidth, 

-  Radar  pulse  repetition  interval  and  pulse  width, 

-  Superhet  time  to  scan  and  time  on  frequency  band, 

-  Mean  hopping  time  and  pulse  width,  etc. 


In  this  thesis  the  above  four  mentioned  window  functions 
will  be  utilized  for  the  analysis  of  the  probability  of 
coincidence. 

Since  the  window  functions  are  periodic  pulse  trains, 
their  coincidences  will  also  be  periodic  [Ref.  2,  p.50J  and 
will  have  the  shape  of  a  pulse  train  also.  The  coincidences 
pulse  train  will  have  a  period  given  by: 


T  = 
o 


where : 


T.  =  Period  of  the  jth.  window  function 
J 

Tj  =  Width  of  the  jth.  window  pulse 


(4-2) 


Computational  formulas  derived  from  the  above  are  found 
in  Appendices  H  and  1  of  this  thesis  for  the  desired  windows. 
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The  mean  duration  of  the  overlaps  (r0)  is  given  by! 


(4-3) 


n  p  1 
)  =  1  l  3 


The  probability  of  coincidence  as  function  of  time,  Pc(t) 
from  Ref.  15  is  expressed  by: 


Pc(t) 


1  -  K  exp 


(4-4) 


where : 


K  =  1  -  Po 


(4-5) 


The  average  coincidence  fraction  Pq  is  given  by: 


n  r  . 

J 


p° '  FT  t 


(4-6) 


The  required  values  for  T^  and  T ^  are  defined  as 


foil ows : 


a)  Antenna  scan  window :  For  the  antenna  scan  window, 

T.  is  the  scan  rate  of  the  antenna  in  seconds  (T  ),  and 
J  S 

r  is  the  time  the  3  dB  beamwidth  takes  to  illuminate  any 
a 

target,  say: 

Ta  =  9  3dE  /  [360°  •  Scan  rate  (°/sec.)  ]  (4-7) 

b.  )  Radar  pulse  train;  The  radar  pulse  train  is  by 
itself  a  window  function  defined  by  its  period  PR1  and 


Pulse  Width  (PW). 
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c.  )  Superhet  scan  and  time  on  frequency :  The  sweep 
period  and  acceptance  bandwidth  of  a  superhet  receiver 
must  be  carefully  chosen  because  they  not  only  affect  the 
probability  of  coincidence  but  also  the  noise  level  of 
the  system.  The  time  the  receiver  dwells  on  the 
frequency  band  of  the  emitter,  r represents  the  window, 
and  to  ensure  interception,  this  time  must  be  greater 
than  the  pulse  repetition  interval. 

Ortiz  and  Giaquinto  [Refs.  9  and  103  have  chosen  T ^  to 
be  5  times  the  radar  PR1  for  calculation  purposes.  The 
time  on  frequency  is  defined  by: 


Ba  +  Bs/2 


D 


Ts  =  5  PR  1 


(4-8) 


where  for  example: 


Ba 

Bs 

D 

Ts 


100  MHz; 
250  KHz; 

2  GHz; 
Scan  time 


Receiver  acceptance  bandwidth 
Radar  frequency  spectrum  app. 
Sweep  bandwidth  of  receiver 
of  receiver 


1/PW 


From  formula  4-7  we  can  calculate  the  values  of: 
Ts  =  0.332  s.  and  r  =  0.05  s;  which  were  used  in  the 
calculations  of  the  Appendices. 
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The  scan  rate  is  also  expressed  as: 


(4-9) 


The  shortest  time  on  frequency  Tf  is  the  minimum 
time  required  by  the  detector  for  energy  buildup,  and  is 
inversely  proportional  to  the  acceptance  bandwidth, 
Tf  (min)  =  1/Ba. 

Therefore  if  this  value  is  replaced  in  equation 
4-19,  we  will  be  able  to  obtain  the  fastest  scan  rate 
achievable  without  loss  of  sensitivity  [Ref.  3], 

D  2 

-  =  Ba  (4-10) 

Ts 
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(4-11 ) 


Ba2  =  1010  MHz/sec 

and  actually  our  sweep  rate  is: 

D  2  GHz 

-  =  -  =  6000  MHz/sec. 

T s  0 • 332  s . 

Which  confirms  that,  for  our  example,  we  are 
scanning  with  a  scan  rate  less  than  the  maximum  allowed. 

d. )  Mean  hopping  time  and  pulse  width:  The  last  window 
function  of  the  model  is  derived  from  the  fact  that  the 
radar  is  hopping  over  five  different  frequencies, 
therefore  according  to  Self  A.  [Ref.  15]  each  frequency 
can  be  treated  as  an  individual  radar.  If  the 
frequencies  are  uniformly  distributed,  the  period  for  a 
particular  frequency  to  be  repeated  on  the  average  is 
equal  to  five  times  the  radar  PRI,  and  since  the  time  on 
frequency  of  the  superhet  receiver  covers  up  to  five 
times  the  PRI,  it  can  intercept  it.  This  problem  must 
be  considered  also  for  an  1FM  receiver  since  for 
correlation  purposes,  the  receiver  needs  to  check  the 
hopping  period.  The  window  width  on  the  average  will 
include  two  radar  pulses  since  the  receiver  is  dwelling 
for  a  time  equal  to  five  times  the  radar  PRI,  and  as 
stated  above,  with  this  particular  kind  of  hopping  the 
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mean  time  for  any  particular  frequency  to  appear  again 
is  five  times  the  PRI. 

C.  PROBABILITY  OF  IDENTIFICATION 

The  identification  of  a  radar  emitter  is  done  by 
comparing  the  received  parameters  with  those  stored  in  an 
internal  library.  Assuming  that  the  stored  parameters  were 
accurately  obtained  and  recorded,  the  probability  of 
identification  is  a  function  of  the  quality  of  the  receiver 
measurements,  the  number  of  parameters  obtained,  the 
characteristics  of  the  processor  and  the  electromagnetic 
density  of  the  environment.  Since  in  this  thesis  we  are 
assuming  a  low  intensity  kind  of  encounter,  we  can  assume  that 
the  processor  will  be  able  to  handle  appropriately  all  the 
information  obtained  by  the  receiver  subsystem.  The  last 
assumption  made  the  identification  problem  a  direct  function 
of  the  Receiver  measurement  accuracy/resolution,  and  the 
number  of  parameters  measured. 

From  information  theory  we  can  state  that  in  order  to 
reduce  uncertainty,  any  decision  must  be  made  over  as  many 
dimensions  as  there  are  available  from  the  information  source. 
This  is  particularly  true  in  EW  since  it  is  required  to  sort 
classify  and  identify  the  signals  intercepted;  and  all  modern 
radar  transmitters  try  to  randomize  their  parameters  in  order 
to  complicate  the  evaluation  process  at  ESM  reception. 


Davies  C.  L.  [Ref.  15,  p.166]  stated  that  the  major 

parameters  required  for  characterizing  an  emission  ares  angle 
of  arrival,  carrier  frequency,  pulse  repetition  interval  (PR1 ) 
and  pulse  width  (PW);  and  the  availability  and  accuracy  of 
them  improves  the  identification  capability  of  the  System. 
From  the  parameters  cited  Angle  of  Arrival  (AOA)  is  the  most 
important  since,  while  the  other  parameters  may  be  varied  by 
the  transmitter,  it  is  impossible  to  mask  its  AOA.  The 
remaining  parameters,  arranged  in  order  of  importance,  are 
frequency,  PRI  and  PW. 

The  parameters  availability  can  be  modeled  by  using 
normally  distributed  accuracy  formulas  corrected  to  obtain 
realistic  performance  that  is  within  0.1  to  1%  of  the  original 
parameter.  Using  one  standard  deviation  (±  1  a)  as 

descriptor  of  variability  gives  a  probability  of  occurrence 
of  only  68.27  X.  If  we  use  two  standard  deviations  (±  2  <x), 
then  the  probability  of  occurrence  raises  to  95.44% 
Accuracy  and  resolution  are  two  closely  related  terms. 
Resolution  is  the  finest  resolving  power  of  the  system  that 
will  allow  for  differentiation  of  two  signals  closely  spaced 
in  the  parameter  spectrum.  On  the  other  hand,  accuracy  of  the 
measurements  depends  on  many  factors  but  its  minimum  value  is 
defined  by  the  resolution  of  the  system.  It  is  important  to 
state  that  a  system  cannot  measure  more  accurately  than  its 
resolving  capability,  but  signal  variations  that  produce 
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inaccuracies  may  well  affect  a  parameter  measurement  if  they 
are  significant;  that  is  bigger  than  the  resolution  cell. 

For  the  case  of  frequency  measurements,  the  superhet  and 
IFM  resolution  problems  are  different.  In  the  first  case,  the 
resolution  is  defined  by  the  IF  bandwidth  of  the  system,  while 
in  most  modern  IFM  systems  it  is  given  by  the  least 
significant  bit  of  the  digital  parameter  word. 

For  the  model  of  this  thesis,  is  assumed  an  10  bit 
frequency  word  in  the  IFM  system  which  makes  the  frequency 
resolution  to  be  approximately  2  MHz.  out  of  a  2  GHz. 
frequency  band. 

The  accuracy  measurement  formulas  derived  from  chapter 
three  of  this  work  and  corrected  for  general  EW  application 
are  as  follows: 

<tF 

Frequency  identification  occurs  with  a  high  probability 

if  (4-o-F)/f  is  <  0.001;  and  if  the  frequency  resolution  is 
c 

also  within  that  limit  (0*001  of  the  original  frequency). 


it  2 -PW3- Bn* SNR 


(4-12) 


0. 35  -  2 

aPRI  = -  (4-13) 

0.8  B  *  SNR'5 
v 


PRi  identification  occurs  if  ( 4 - c - PR  I ) /PR  1  '  0.01 
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<xPW 


(4-14) 


.  7 

B  -SNR'5 
v 

PW  identification  occurs  if  (4-a-PW)/PW  <  0.01 

The  probability  of  identification  can  be  modeled  from 
combinatorial  analysis  and  using  the  utility  factors  approach 
recommended  by  Glenn  [Ref.  73.  Therefore  we  will  be  able  to 
weight  properly  each  of  the  possible  combinations  of  the 
parameters  at  reception  according  to  the  importance  of  each 
of  them  and  make  use  of  MATHCAD’ s  MAX  function  to  select  the 
most  valuable  combination.  The  utility  factors  of  the 
parameters  at  reception  are: 


Angle 

of  arrival 

wl  =  0.3 

Carrier  frequency 

w2  =  0.28 

Pu  1  se 

repetition  interval 

w3  =  0.25 

Pulse 

width 

w4  =  0.19 

A  sensitivity  analysis  was  applied  to  these  utility 
factors  in  order  to  check  their  influence  in  the 
identification  process.  This  will  be  discussed  in  the  next 
chapter  during  the  analysis  of  the  output. 

When  a  parameter  is  measured,  and  its  deviation  is  within 
the  required  limits,  then  we  will  be  able  to  use  it. 
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This  is  modeled  by  using  MATHCAD's  Heaviside’s  step 
function : 

d(x)  =  5  (limit  -  (4<r/p))  (4-15) 

Therefore  if  the  ratio,  deviation/mean,  is  less  than  the 
required  limit  a  "1"  is  returned,  otherwise  it  will  return  a 
"0”.  Finally,  with  the  use  of  combinatorial  analysis,  an 
identification  array  is  modeled  that  includes  all  the 
available  parameter  combinations  at  reception  that  may 
identify  a  certain  transmitter,  weighting  each  particular 
combination  according  to  which  parameters  it  is  composed  of, 
and  using  the  standard  deviation  model  to  determine  which 
parameters  were  available  at  every  measurement.  In  this  way, 
the  combinations  are  arranged  from  top  to  bottom  with  the  most 
valuable  at  the  top  and  the  least  at  the  bottom.  The 

weighting  factor  used  in  every  row  (combination),  was  chosen 

to  be  the  square  root  of  the  summation  of  the  weights  of  the 

respective  parameters  involved.  The  usage  of  this  array 
guaranteed  that  only  the  row  which  included  all  the  available 
measurements  will  be  selected,  avoiding  redundant  choices. 

This  is  true  because,  for  example,  if  three  parameters  are 
identified,  ail  combinations  will  acquire  a  certain  value  but 
only  the  one  which  includes  all  three  must  have  been  selected. 
Initially  the  weighting  factors  were  defined  just  by  the  sum 
of  the  weights  of  the  respective  deviations  within  each 
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combination;  but  this  approach  was  discovered  to  overweigh  the 
combinations  which  had  more  elements  (or  more  weight). 
Therefore  when  only  two  or  three  parameters  were  identified 
the  probability  of  identification  was  decided  not  by  the 
combination  which  contained  the  identified  parameters  but  by 
the  more  valuable  combination  (the  one  which  includes 
measurements  of  all  the  parameters)  even  though  it  did  not 
have  all  the  deviations  identified.  This  was  because  the 
weights  overcame  the  lack  of  parameters.  It  is  clear  then 
that  the  model  was  not  correct,  since  even  though  it  included 
all  the  variables  of  the  problem,  it  always  selected  the  same 
combination  (the  one  with  the  highest  weighting  factor). 

(wl+w2+w3+w4)5  -0.25  ( dAOA+dF+dPR I +dPW ) 

( wl + w2+w3 )  ’5  -0.333  ( dAOA+dF+dPR 1 ) 

( wl +w2+w4 )  ‘5  .0.333  ( dAOA  +  dF+dPW ) 

(wl+w3+w4),S  .0.333  ( dAOA+dPR I +dPW ) 

(w2+w3+w4)'5  .0.333  (dF+dPRI +dPW) 

(wl  +  w2)'5  -0.5  (dAOA  +  dF) 

(wl+w3)'5  -0.5  (dAOA  +  dPRI  )  (4-16) 

(wl+w4)'5  *0.5  (  d  AOA  +  dPW ) 

( w2  +  w3 ) ’ 5  *0. 5  (dF  +  dPRI) 

(w2  +  w4)*5  .0.  5  (dF  +  dPW) 

(w3  +  w4)‘5  *0.  5  (  dPRI  +dPW) 

(w2)'5  •  dF 
(w3)'5  •  dPR I 
(w4)’5  •  dPW 
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The  weight  differences  from  row  to  row  were  also  too 
significant.  As  a  result,  weight  difference  became  the  key 
in  the  selection  process  even  though  in  some  cases  there  was 
a  lack  of  other  information. 

In  order  to  solve  this  problem,  the  differences  between 
the  rows  of  the  array  created  in  the  original  model  were 
reduced  by  taking  their  square  root.  Then  by  using  MATHCAD’ s 
MAX  function,  we  were  able  to  chose  the  combination  which  had 
the  highest  value  in  the  probability  of  identification  array. 

The  following  example  illustrates  the  weighting  factors 
problem.  Suppose  that  the  system  has  identified  frequency, 
PRI,  and  PW,  and  MATHCAD’ S  MAX  function  is  used.  Then: 

dAOA  =0;  dF  =1;  dPRI  =1;  dPW  =  1 

a)  Method  No .  1_:  Weight  =  Sum  of  weights 

1.0  (0.25)  ( dAOA+dF+dPR I +dPW)  =  .75  .  .  wrong  selection 

0.9  (0.33)  ( dAOA+dF+dPR 1 )  =  .594 

0.7  (0.33)  ( dAOA+dPR I +  dPW )  =  .462 

0.6  (0.33)  (dF+dPRI +dPW)  =  0.6 

.  5 

b)  Method  No . 2 :  Weight  =  (Sum  of  weights)' 

1.0  (0.25)  (dAOA+dF+dPRl+dPW)  =  .75 

0.948  (0.33)  (dAOA+dF+dPR I )  =  .626 

0.637  (0.33)  ( dAOA+dPR I +dPW )  =  .557 

0.774  (0.33)  (dF+dPRI +dPW)  =  .774  .  .  right  selection 
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Therefore  by  using  MATHCAD’s  MAX  function  in  method  No. 2 


the  correct  selection  was  made. 


Parameters  Measured 

PROBABILITY  OF 

AOA 

FREQ 

PR  1 

PW 

IDENTIFICATION  P(id) 

X 

X 

X 

X 

1.00 

X 

X 

X 

0.911 

X 

X 

D 

0.877 

X 

X 

X 

0.86 

X 

X 

X 

0.35 

X 

X 

0.761 

X 

X 

0.741 

X 

X 

.7 

X 

X 

0.728 

X 

X 

0.685 

X 

X 

0.  663 

X 

0.529 

X 

0.5 

X 

0.  436 

Table  1',  Probability  of  Identification  based  on  combinatorial 
analysis,  weighting  factors,  and  parameter  accuracies. 

The  combination  that  matched  the  parameters  identified 
was  weighted  the  most  valuable  compared  with  the  other 
possible  combinations  within  the  array. 
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The  angle  of  arrival  identification  was  particularly 
difficult  to  model  because  it  depends  on  the  direction  finding 
(PF)  receiver  system  employed.  Consequently,  in  order  to 
solve  it  we  established  an  AOA  rms.  error  of  4  degrees  as  a 
reference  and  -50  dBm  as  its  sensitivity  level.  These  are 
typical  values  of  performance  by  equipment  used  for  tactical 
purposes.  Generally,  most  ESM  systems  have  a  separate  DF 
receiving  system  that  is  interfaced  with  the  radio  frequency 
system.  The  most  common  DF  system  uses  amplitude  comparison 
techniques  and  its  receivers  are  broadband  crystal  video 
in  an  omni-directional  array  of  4,  6  or  8  antennas. 

The  DF  accuracy  obtained  in  the  model  will  be  a  function 
then  of  the  strength  of  the  signal  at  reception,  the 
sensitivity  level,  and  the  ratio  between  the  actual  angular 
error  of  the  system  compared  with  the  stated  reference.  The 
actual  rms  error  was  set  to  5  degrees. 

As  documented  in  this  chapter,  the  POI  is  a  function  of 
three  independent  factors,  probability  of  detection, 
probability  of  coincidence,  and  probability  of  identification. 
Each  one  of  these  factors  being  particularly  defined,  the 
measure  of  effectiveness  reflects  the  complex  interaction  of 
the  many  inputs  involved,  and  allows  the  designer  to 
selectively  analyze  the  parameters  influence  in  the  overall 
performance  variation. 
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V.  ANALYSIS 


In  the  analysis  section  of  this  thesis  we  will  go  over 
the  outputs  of  the  different  MATHCAD  files  which  were 
developed  in  order  to  represent  the  P01  as  explained  in 
chapter  four.  The  causes  and  effects  related  to  each 
particular  set  of  curves  will  be  discussed.  The  generated 
curves  are: 

-  Main  lobe  and  side  lobes  probability  of  detection 

-  Main  lobe  and  side  lobes  probability  of  coincidence 

-  Main  lobe  and  side  lobes  probability  of  identification 

-  Main  lobe  and  side  lobes  probability  of  intercept 

The  curves  were  modeled  from  both  IFM  and  superhet 
receiver  special  characteristics,  as  defined  in  appendices  B 
and  C.  The  aircraft  was  initially  placed  at  4000  m.  height, 
(279  Km.  maximum  range)  on  the  ESM  horizon  line  (tangential 
to  the  Earth),  beyond  the  estimated  radar  maximum  range,  which 
will  allow  us  to  observe  typical  performances. 

A.  PROBABILITY  OF  DETECTION  ANALYSIS 

The  analysis  starts  with  the  probability  of  detection 
curves  Figures  S,  10  and  11.  Figure  9  shows  the  probability 
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of  detection  of  the  radar's  sidelobes  as  a  function  of  tame 
for  a  superhet  receiver.  Each  time  unit  represents  200 
seconds  of  flight.  Therefore  after  flying  000  sec.  (app.  at 
200  Km.  from  the  target),  this  receiver  will  act  as  an 
omnidirectional  antenna  by  receiving  the  radar  all  the  time. 
It  should  be  noted  that  the  probability  of  detection  for  the 
radar  main  lobe,  as  expected,  was  always  unity  across  the  time 
window  of  interest.  The  equations  used  from  Tsui  [Ref.  3,  Ch. 
2],  which  include  the  processing  gain  effect  due  to  the  ratio 
Br/Bv,  actually  reduce  the  SNR  required  at  the  receiver  to 
obtain  a  certain  probability  of  detection  given  a  probability 
of  false  alarm.  In  fact,  the  required  SNR  to  produce  a  0.999 
probability  of  detection  is  approximately  10  dB.,  when  the 
Pfa  =  that  might  seem  too  small. 


0  tl  2200  seconds 

t 

279  Km.  62  Km. 

» 

Fig.  9  Radar  side  lobes  probabiol i t 1 y  of  detection  vs.  time 
superhet  receiver  case. 
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The  probability  of  detection  versus  time,  for  the  example 
1FM  receiver,  showed  that:  due  to  the  rms.  noise  power 
increase  (wider  bandwidth),  this  receiver  was  unable  to  detect 
the  radar  sidelobes  during  all  the  flight  time.  Near  the  end 
of  the  time  line  (62  Km.  from  target)  the  SNR  for  the  radar 
side  lobes  was  just  0.88  dB;  too  low  for  achieving  detection, 
and  this  agrees  with  the  curves  from  Tsui  tRef.  3,  Ch.  23. 
In  any  case,  as  Figure  11  shows,  it  did  detect  the  radar 
mainlobe  with  a  probability  of  one  along  the  trajectory. 


Fig.  10  Radar  main  lobe  probability  of  detection  vs.  time; 
superhet  receiver  case. 

The  analysis  of  the  probability  of  detection  factors 
showed  that  the  superhet  receiver  had  the  advantage  of  having 
detection  of  both  the  radar  main  and  side  lobes.  On  the  other 
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hand  the  IFI1  receiver  was  able  to  detect  the  radar  emission 
only  when  the  radar  main  lobe  pointed  at  it. 


PdmIFM 


279  Km. 


2200 


62  Km. 


Fig.  11  Radar  main  lobe  probability  of  detection  vs.  time; 
I FM  receiver  case. 


B.  PROBABILITY  OF  COINCIDENCE  ANALYSIS 

The  probability  of  coincidence  vs.  time  in  the  superhet 
receiver  was  very  low  as  can  be  obserrved  in  Figure  12,  where 
after  2200  seconds  of  flight  the  probability  of  coincidence 
is  0.06.  This  is  the  disadvantage  of  using  a  narrowband 
scanning  receiver  against  a  frequency  hopping  radar. 

From  the  respective  I1ATHCAD  file,  the  average  time 
between  coincidences  (period)  from  a  superhet  receiver  for  the 
mainlobe  and  sidelobes  was  found  to  be  8  h.  50  min.  30  sec., 
and  3  min.  42  sec.  respectively  (Tom  and  1  os  in  Appendix  H). 


Fig.  12  Superhet  receiver  main  lobe  probability  of 
coincidence  versus  time. 


Fig.  13  Superhet  receiver  side  lobes  probability  of 
coincidence  versus  time. 
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The  low  performance  induced  by  increasing  the  number  of 
window  functions  is  drastic,  which  explains  why  modern  radar 
designs  are  trying  to  randomize  (within  limits)  most  of  the 
transmitted  parameters. 

One  way  to  improve  the  probability  of  coincidence  is  by 
elimination  of  one  window  function  or  optimization  of  the 
window  functions  that  are  under  our  control;  the  first 
solution  is  possible  if  1FM  receivers  are  used.  IFM  receivers 
can  observe  instantaneously,  wider  bandwidths  eliminating  the 
receiver-scan  window  function  (setting  it  equal  to  one).  This 
effect  is  shown  in  Figures  14  and  15,  for  main  and  side  lobes 
probability  of  coincidence,  respectively. 


Fig.  14  IFM  receiver  main  lobe  probability  of  coincidence 
versus  time. 
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In  Figure  15  a  0.9  probability  of  coincidence  (mainlobe) 
at  the  end  of  the  time  window  of  interest  is  achieved;  and  in 

Figure  16  it  can  be  seen  that  the  probability  of  coincidence 

approaches  unity  very  fast,  which  is  a  tremendous  improvement 

compared  with  the  superhet  case. 


1 


Pcs  1 1 FM 

t 


0 

0  tl  40 

t 

279  Km  275  Km. 


»* 

1* 

_  __ 

»* 

f 

1 

Fig.  15  IFM  receiver  side  lobes  probability  of  coincidence 
versus  time. 

In  general  it  can  be  said  that  the  probability  of 
coincidence  curves  gave  the  IFM  receiver  a  significant 
preformance  advantage  when  compared  with  the  superhet  receiver 
performance  against  both  the  radar  main  and  side  lobes. 
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C.  PROBABILITY  OF  IDENTIFICATION  ANALYSIS 

The  probability  of  identification  is  quite  empirical 
since,  even  though  the  accuracy  of  the  measurements  relates 
to  identification  of  parameters,  these  formulas  were  evaluated 
with  assumed  errors  which  try  to  mimic  realistic  ESM 
scenarios. 

For  the  case  of  superhet  receiver,  it  can  be  observed  in 
Figure  16  that  for  main  lobe  reception,  the  probability  of 
identification  has  a  value  of  0.78  during  all  the  time  of 
interest.  This  is  explained  because  the  model  is  identifying 
all  the  parameters  except  frequency  (wider  resolution). 
Additionally  AOA,  which  was  set  to  be  the  ratio  of  a  reference 
value  to  the  actual  accuracy,  had  an  accuracy  of  0.8  (4/5 
=  0.8);  therefore,  the  identification  array  chose  the  fourth 
combination  giving  us  a  probability  of  identification  of  0.79. 

For  the  superhet  side  lobes  case,  the  probability  of 
identification  has  a  value  of  0.48  from  t  =  0  s.  to  1450  s. 
where  it  jumps  to  0.663.  The  first  value  is  due  to  PR  1 
measurement  alone,  and  after  closing  in  range  for  1450  s.  the 
system  is  able  to  identify  additionally  the  emitter  PW,  thus 
improving  its  identification  contribution. 

For  an  I FM  receiver,  the  mainlobe  identification  problem 
is  depicted  in  Figure  17,  where  (at  t  =  1320  s.)  the 
probability  of  identification  jumps  from  0.84  to  0.95.  This 
jump  is  due  to  PW  accuracy  improvement  at  that  particular 


57 


time,  providing  thereafter  all  the  parameters  needed  for 
identification  purposes. 


identification  versus  time. 


Fig.  17  Superhet  receiver  sidelobes  probability  of 
identification  versus  time. 
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279  Km.  62  Km. 


Fig.  18  1FM  Main  lobe  probability  of  identification  vs,  time 

Figure  19  plots  the  side  lobes  case  of  the  same  problem. 
In  here  we  find  a  straight  line,  representing  a  constant 
probability  of  Identification  value  of  0.728.  This  value 
comes  from  having  achieved  only  PR1  and  frequency  accuracy, 
no  other  parameter  was  available  during  the  time  window. 

The  probability  of  identification  plots  have  shown  an 
considerable  advantage  of  the  IFM  on  main  and  side  lobes 
emitter  identification.  The  main  source  of  this  difference 
is  due  to  the  frequency  identification  capability  of  the  IFM 
receiver  (fine  frequency  resolution).  Frequency  availability 
has  been  weighted  as  one  of  the  most  valuable  parameters  for 
identification  purposes. 
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Fig.  19  1FM  side  lobes  prob.  of  identification  vs.  time. 

D.  POI  ANALYSIS 

Finally  all  the  above  probabilities  are  used  to  model  the 
overall  probability  of  intercept  (POI),  this  may  be  expressed 
as  A  function  of  two  mutually  exclusive  events: 

POI  =  POIml  +  POIsl  -  (POIml -POIsl )  (5-1) 

where : 

POIml  =  Probability  of  intercept  main  lobe 
POIsl  =  Probability  of  intercept  side  lobes 

Each  one  of  the  above  variables  are  derived  from  three 
mutually  independent  events: 
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POlml 

zz 

PDm 1  PI  Dm!  PCml 

(5-2) 

POI  s  1 

= 

PDs 1  PIDsl  PCs  1 

(5-3) 

where  : 

PDm  1 

= 

Probab i 1 i ty 

of 

detection  -  main 

1  obe  - 

PDs  1 

= 

Probab i 1 i ty 

of 

detection  -  side 

lobes  - 

P  I  Dm  1 

= 

Probab i 1 i ty 

of 

identification  - 

main  lobe 

- 

PIDsl 

= 

Probab i 1 i ty 

of 

identification  - 

side  lobes 

- 

PC-m  1 

- 

Probab i 1 i ty 

of 

coincidence  -  main  lobe  - 

PCs  1 

= 

Probab i 1 i ty 

of 

coincidence  -  side  lobes  - 

The  above  data  was  generated  for  every  time  step  within 
each  of  MATHCAD  file.  The  outputs  were  also  stored  in  the 
form  of  files  in  such  a  way  that  the  P01  file  recovered  the 
component  files,  did  the  calculations  above  explained,  and 
outputed  the  POi  results.  The  plot  shown  in  Figure  19 

represents  the  overall  probability  of  intercept  for  both  IFM 
and  superhet  receivers  in  this  particular  scenario. 

The  most  significant  aspect  about  Figure  20  is  that  the 
IFM  performance  did  not  make  use  of  any  Side  lobe  interception 
since  it  was  never  able  to  detect  it.  In  addition,  the 
performance  of  the  superhet  receiver  was  not  as  good  as  that 
of  the  IFM  even  though  it  achieved  very  good  detection  (SNR) 
and  identification  (theoretical  from  SNR)  of  main  and  side 
lobes  along  the  time  line.  From  this  we  can  deduce  that 
the  probability  of  coincidence  is  the  most  significant  factor 
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Fig.  20  Probability  of  intercept  versus  time,  superhet  and 
1FM  receivers  performance. 


in  these  scenarios  and,  as  has  been  discussed  in  many 
publications,  it  represents  a  key  point  for  superhet  receivers 
which  suffer  the  larger  degradation  in  their  interception 
capabilities.  In  order  to  investigate  this  further,  it  was 
decided  to  test  the  model  against  the  same  radar  but  while 
transmitting  a  constant  carrier  frequency.  The  probability 
of  intercept  plot  obtained  is  in  Figure  21. 

This  confirms  the  influence  of  the  probability  of 
coincidence  in  the  general  performance.  By  fixing  the  carrier 
frequency,  one  eliminates  one  of  the  window  functions,  thereby 
improving  both  systems  in  an  equal  manner.  Now  the  superhet 
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Fig.  21  Probability  of  intercept  vs.  time,  superhet  and  IFM 
performances  against  a  fixed  frequency  radar. 

receiver  performs  better  until  app.  t  =  1300  s.  because  the 
additional  weighting  of  having  both  main  and  side  lobes 
interception  supplies  a  significant  increase  in  POI.  This 
overcomes  the  additional  window  function  that  the  Superhet  was 
encountering  (scanning  in  frequency). 

Both  curves  show  a  "jump"  approximately  between  t  =  1300 
and  1450  s.  This  is  derived  from  the  PW  identification,  as 
explained  above. 

In  general  terms,  when  dealing  against  a  fixed  frequency 
radar,  the  superhet  receiver  experiments  a  considerable 
improvement.  This  improvement  is  basically  due  to  better 
probability  of  coincidence  performance,  and  if  frequency 
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resolution  is  refined  by  the  use  of  a  narrower  band,  better 


identification  values  can  be  obtained. 

E.  SENSITIVE  ANALYSIS  OF  WEIGHTING  FACTORS 

Finally,  the  weighting  factors  were  tested  in  order  to 
check  their  influence  in  the  probability  of  identification  and 
in  the  probability  of  intercept.  The  weighting  factors  as 
summarized  in  Table  2.  show  the  use  of  the  following  weights: 


Figure  21 

Weight  set  a. 

Figure  22 

Weight  set  b. 

Figure  23 

Weight  set  c. 

dAOA  =0.4 

0.35 

0.35 

dF  =  0.3 

0.35 

0.  3 

dPRI  =  0.2 

0,  25 

0.  2 

dPW  =  0. 1 

0.05 

0.  15 

Table  2.  Weighting  factors  sets. 


As  can  be  seen  in  Figure  21,  the  use  of  the  first  set  of 
weights  does  not  affect  the  initial  probability  of  intercept 
(Figure  19)  except  by  a  small  reduction  in  the  superhet 
receiver.  The  performance  is  otherwise  nearly  the  same  as  for 
Figure  20.  Reviewing  the  probability  of  identification  plots 
it  was  found  that  the  side  lobes  jump  at  t  =  1450  s.  changed 
initially  from  0.479  to  0.648  (Figure  16)  but  in  Figure  21 
they  changed  from  0.447  to  0.557,  For  the  IFM  case  the 
variations  were  not  significant. 
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Fig.  22  Superhet  and  IFM  performances;  dAOA  =  0.4,  dF  =  0.3, 
dPR I  *  0.2,  and  dPW  =0.1 

Following  the  a-set  weights  evaluation,  b-set  weights 
were  substituted  and  the  plot  of  Figure  22  was  obtained.  It 
can  be  noted  a  significant  decrease  in  the  probability  of 
intercept  curves,  and  especially  in  the  superhet  curves,  as 
compared  with  Figure  20;  but  if  they  are  compared  with  Figure 
22  it  can  be  found  after  the  jump  both  Figures  show  no 
difference.  One  observes  an  approximate  0.2  improvement,  in 
the  superhet  receiver,  before  the  jump  takes  place.  With  the 
last  weights  set,  Figure  23  is  obtained;  which  is  very  similar 
to  Figure  22. 


65 


The  only  difference  is  that  the  jump  in  the  superhet  POI 
is  larger  with  set  c.  This  is  because  PW  is  now  weighted  more 
heavi 1 y. 

The  precedence  used  in  this  work  for  the  parameters 
weighting  should  be  kept  unless  extremely  good  processor 
performances  would  allow  the  use  of  higher  utility  factors 
consequently  modifying  the  established  precedence.  It  would 
be  desirable  also  to  test  operationally  fcr  the  values  of  each 
parameter  with  special  consideration  of  the  processor 
capabi 1 i t i es . 

As  can  be  observed  in  the  plots,  the  shape  of  the  POI 
curves  would  not  be  altered  but  their  time  to  achieve  a  POI 
value  could  be  reduced  with  variations  in  the  weighting 
factors.  In  any  case,  the  approach  used  in  this  thesis  is 
useful  since  it  allows  good  approximate  results. 
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VI.  CONCLUSIONS 


A.  GENERAL  CONCLUSIONS 

1)  During  ESM  systems,  performance  evaluation  many 
variables  and  complex  trade  offs  between  them  must  be 
considered.  The  model  used  in  this  work  is  robust  and  expands 
the  qualitative  evaluations  commonly  done  with  a  "good", 
a  ,  r  '  •  o  '  cor  f  1  v  -  a  d  i  n  g  p-o .  i  c  /  » 

Z  Within  the  scope  and  scenarios  of  this  thesis  it  can 
oe  concludes  from  the  measure  of  effectiveness,  that  the  IFM 
receiver  had  a  better  performance  thatn  the  Superhet. 

.  v  it.  *  Ota  .  iMu  w  1 1  i *  ^  e.j  li  o r  c  /  ag  i  i  e  radar-.  the  fact 
r<f  ha  v  i  r  r.  a  ve*-/  sensitive  receiver  by  itself  does  not  solve 
if.f.  ;  .  te  r  rr.px  ic.f  problem.  as  indicated  by  the  Superhet 
rece.ve-  cu>  ve  .  The  ESN  designer  must  work  within  the 
coi.'.;  i-.cri.e  c-rcolem  to  imprive  tne  probability  of  coincioence, 
modi J  ■  j  ng  ire  wi'TOvs  uno6r'  his  control.  Rotating  antennas 
a  *=  .  ~t  f?rroiir;f*.s  *  o«*  trcL-  cel  purposes  since  toeir  use  will 
imply  an  e<tru,  window  function  in  the  coincidence  problem. 


B.  SPECIFIC  CONCLUSIONS 


1.  Tne  per  f  crmao.-.e  oose-  -ei  in  the  super  net  receiver 
aoev.nst  a  f  i  >  e  --  frequency  transmitter  was  superior  to  the  IFM 
c?  >“i  '•  -t  r  ■  defm-t  r.-  tne  sioelnnes.  I 


v  i 


r  *  detm-t 


t  is 


well  1  r.owr,  that  superhet  receivers  are  capable  of  achieving 
t<Isc  hignef  frequency  accuracy  as  compared  with  IFM  receivers. 

2>  The  probability  of  coincidence  requirement  for  a 

tactical  ESh  system  must  be  always  considered. 

3)  Sensitivity  still  remains  a  very  important  parameter 
since  signal  value  is  critical  at  every  step  in  the 
ca 1 cu 1 ations . 

C .  OPERAT I  ONAl  CONCLU‘3 1 ONS 

1)  In  peace  time.  when  threat  radars  are  more 

iriGiscrete,  a  superhet  receiver  is  the  optimum,  choice. 

2  In  wartime  tactical  operations,  the  IFM  is  superior 
because  cf  its  shorter  time  to  interception  even  though  its 
accuracies  are  not  as  high.  For  initial  tactical  success,  the 
£:■!'  opt  ■  a \ci  •*«.  st  utilace  ns.  cssets  tc  the  maximum  continuous 
e-  tert  pc-ssiols. 

3  t  Df'  i  tu  tac t i ca 1  surveillance  operations  the  aircraft 

snc-uld  Siwfvf  f  1  be 1  ow  the  line  of  sight  of  tne  emitter  which 
is  tangential  tbs  earth's  surface  sno  awa-  from,  the 

langencsal  poire.  Tnis  is  the  worst.  Detection  region  for  the 
radar,  and  if  the  aircraft  is  flying  high  enough,  the  E5M 
range  oecom.es  tne  best  weapon  at  hand.  As  the  fcSM  aircraft 
dc.eu  its  approach  towards  the  emitter,  it  must  .ower  its 
oititude  continuous!  -  ;n  order  to  stay  ir.  the  radar 
diffractic  rone.  ES r^roer  cbtainable  with  sensitivities 


of  -65 

maxim  tiff. 


■  dBm.  or  better  are  considerably  greater  than  the 
radar  detection  range. 


VII.  RECOMENDAT IONS 


i;  Operational  test  should  be  conducted  to  compare  and 
validatethe  data  presented  here.  The  combinatorial-weighting 
factors  approach  presented  appeared  to  give  artificially  high 
accuracy  values. 

It  During  peace  time  operations;  ESM  aircraft  crews 
should  concentrate  their  efforts  toward  fine  tuning  radar 
signals  in  oroer  to  upgrade  current  threat  libraries. 


APPENDIX  A:  SENSITIVITY  REQUIREMENTS 


SENSITIVITY  REQUIREMENTS  FOR  AN  ESM  RECEIVER 


The  required  sensitivity  is  a  function  of  the  enemy  radar  and 
of  the  operational  requirements  imposed  such  as  detection  range 
or  pulse  density. 


Enemy  radar  parameters: 
3 

P  :=  250  •  10  watts 

t 

G  :=  1585 
t 

-6 

t  :»  4-10  s 
X  :  =  0. 1  m 


L  :  =  2 
t 

g  :=  1 

h  :=  20 
t 


Peak  power 

Antenna  Gain  (app.  32  db) 

Pulse  width 

Wavelength  (in  meters) 

S  band  radar 

Transmitter  losses 

Multipath  coefficient 
Radar  antenna  height  (m) 


ESM  parameters: 

G  :  =  1 

esm 

0  : =  20 

L  :  =  2 

esm 

h  : =  1000 

esm 


Omnidirectional  Antenna 

Aircraf  t’ s  R.  C.  S. 

ESM  receiver  losses 

ESM  aircraft  height  (m) 
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1.-  ESM  AND  RADAR  SENSITIVITY  REQUIREMENTS  AS  FUNCTION  OF  RANGE 


R  :=  1,5  ..1000  Range  variations 

Signal  power  at  Radar  receiver  : 

2  2  4 

P  -G  *  o’  *  X  *  g 
t  t 

S  CR)  :=  10  log  - 

t  3  4 

(4 '  it)  -(1000-R)  -L 

t 

2  2  4 

P  G  cr-X  g 
t  t 

S  (R)  j=  -  (in  Watts) 

tl  3  4- 

( 4  '  tt  )  (1000-R)  -L 

t 


|  (in  dBw) 


dBw 

-110 


S  (R) 
t 


Power  at  Radar  receiver  vs.  Range 
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ESM  noise  floor  app.  -98  dBw  for  an  I FM 
Radar  noise  floor  app.  -125  dBw  (typical) 


RADAR  AND  ESM  MAXIMUM  DET.  RANGES  AS  FUNCTION  OF  SENSITIVITY 


a.-  Radar  maximum  detection  range 


k  :=  1.3810 


J/deg 


Boltzman’s  constant 


To  :=  290 


Kelvin 


B  :  =  - 


Receiver  Bandwidth 


SNR  :=  10,10.5 

min 


.35.5 


Minimum  SNR  reqd.  for 
detection  (10  to  15  dB) 


F  :=  4 


Noise  Figure  (6  dB ) 


i  [SNR  *1  :• 

rmax  min  J 


2  2  4 

P  G  '  X  ’  >7  ■  g 
t  t 


(4  •  tt)  '  k  ’ 


To-B  - F ' SNR  | 


!  i*S  NR  1 

rmax  *u  min  J 


35.5 


10  dB 


15.5  dB 


Radar  maximum  detection  range  as  function  of  SNR 
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b.-  ESM  receiver  maximum  detection  range  (this  is  minimum 
signal  required  above  the  noise  floor  to  in  order  to  generate  a 
detection. 


Sensitivity  range  -105  to  -65  dBw  (-75  to  -35  dBm): 

-11  -10  -7 


S  :=  3.162-10  ,2010  ..3.162-10 

min 


-11  S  -7 


3. 162  10  min  3.14910  watts 

(-105  dB?  -75  dBm)  (-65  dB ;  -35  dBm) 


ESM  performance  in  range  for  detection  of  a  radar  main  lobe  as 
function  of  various  sensitivity  levels. 


2.-  FLIGHT  HEIGHT  (ESM  Rmax  =  Radar  Horizon) 


h 


t 


Radar  height  (Km) 


h 

esi 

h  -  E  S  M  aircraft  heigth  (Km) 

esmKm  1000 


Sensitivity  range,  -iOS  dBw  (-75  dBm)  to  -75  dBw  (-45  dBm): 


H 

1 

H 

H 

1 

- 

s 

o 

— i 

O) 

o 5 

o 

H 

(D 

H 

CO 

tl 

..3.16-10 

min 

Radar 

Horizon  (Km)  equation: 

r  .5 

r  -*5’i 

R 

:=  130.34  ■  j  fh  1  + 

fh  i 

Hz 

LL  esmKm  J 

L  tKm  J  J 

then  the  height  in  Kilometers  will  be  found  from: 


h  fs 

flgt  l  min 


7.672-10 


h 

■■  tKm 


h 

flgt 


Km. 

1000 


10 


-11  S 

3.1610  min 

(-105  dBw;  -75  dBm) 


-8 

1.010  watts 
(-80  dBw;  -50  dBm) 


Flight  height  (in  Km)  for  detecting  the  Radar’s  main  lobe 
(Gt  =  1585)  vs.  various  Sensitivity  levels. 


For  side  lobes  detection,  (Gt  =  0.5) 


-12  -12  -10 
S  :=  1-10  ,210  ..3.16227-10 

min 


h 

f  lgt2 


r» .  1 

L  min  J 


7.672-10 


h 

f  lgt2 


Km. 

10 


.  2  ' - 

-12 

S 

-10 

10 

min 

3.1610  watts 

(-120  dB; 

-90  dBm) 

(-95  dB;  -65  dBm) 

Flight  height  (in  Km.)  for  detecting  the  Radar’s  sidelobes 
(Gt  =  0.5)  vs.  various  sensitivity  levels 


APPENDIX  B:  SETUP  (SUPERHET) 


Superhet  Receiver  case 

In  this  file  the  user  must  state  the  required  parameters  for 
both  the  ESM  receiver  and  the  enemy  Radar  to  be  encountered. 

The  follow  on  calculations  will  automatically  update  their 
output  every  time  a  change  is  made. 


SH1PB0RNE  ENEMY  RADAR  (AIR  SURVEILLANCE) 


hr  :=  20 
Pt  :=  250000 
9 

ft  :=  3-10 

8 

3  •  10 

Xt  - 

ft 

Gtm 1  :=  1585 
Gts 1  : =  .5 
PRF  :=  300 
1 

PR  I  :=  - 

PRF 

-6 

PW  : =  U  ■  10 
Lt  :=  3 
g  :=  1 
RPM  :=  6 
83dB  : =  2.5 
60 

Ta  :=  - 

RPM 


Enemy  radar  height  (m) 

Peak  power  (watts) 

Center  frequency  (Hz) 

Transmitter  wavelength 

Gain  Radar  Main  Lobe  (  32  dB) 
Gain  Radar  Side  lobe  (  -3  dB) 
Pulse  Repetition  Frequency  (Hz) 

Pulse  Repetition  Interval 

Pulse  width  (sec) 

Radar  transmission  losses 
Multipath  factor  0  <  g~2  <  4 
Antenna  Scan  Rate 
Antenna  3dB  Beamwidth 

Period  of  Antenna  Scan 
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63d  B 

•ra  s  - - Ta 

360 


ESM  RECEIVER  PARAMETERS 


Vmph  :=  220 
Vm  : =  Vmph • . 447 
hesm  :=  4000 

9 

FL  :=  2.0  10 

9 

FU  4.0-10 
D  :=  FU  -  FL 
6 

Br  :=  50-10 

6 

Bv  :=  10  10 

6 

Ba  :=  50-10 
Br 

r  :  =  — 

Bv 

gamma  : =  .8 

gamma 


Fn  :=  10 

Fn 

10 

FN  :=  10 
Gesm  :=  1 

-11 

Pfa  :=  l  10 


Time  Antenna’s  Main  Lobe  Points 
in  a  particular  direction 


ESM  platform  speed  (knots) 

ESM  platform  speed  (m/seo) 
Height  of  ESM  platform  (m) 

Minimum  Frequency  (Hz) 

Upper  Frequency  (Hz) 

Frequency  Band  of  Coverage 

RF  Bandwidth  (Hz) 

Video  Bandwidth  (Hz) 

Acceptance  Bandwidth  (Hz) 

Ratio  of  Br  to  Bv 

Integration  efficiency  coeff. 
Processing  Gain 

Receiver  Noise  Figure  (in  dE) 

Receiver  Noise  Figure 
ESM  Antenna  Gain 
Probability  of  False  Alarm 
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Learn  t=  3 


Receiver  internal  losses 


Ts  i  = 


5  D 


•  PRF 


•rf 


1 

Ba  +  — 
PW 


D 


•Ts 


d 


PW 

2 


Time  to  scan  Frequency  band 
(sec.  ) 


Time  a  particular  frequency 
remains  in  the  Receiver  Passband 


Minimum  coincidence  duration 


T  :=  10 


Time  increment  (sec) 


tf  :=  220 


Increment  limit 


Earn  and  Radar  input  parameters  are  written  into  output  files 
called  ESM  and  RADAR  respectively. 


i  :=  1  . . 14  j  :=  1  . .22 


Radar  : = 


i 


hr 

pt 

ft 

xt 

Gtm  1 

Gts  1 

PRF 

PR  I 

PW 

z 

83dB 

Ta 

■ra 

l.t 

ESM  :  = 

j 


Vmph 

Vm 

hesm 

FL 

FU 

D 

Br 

Bv 

Ba 

r 

gamma 

PG 

Fn 

FN 

Gesm 

Pf  a 

Lesm 

Ts 

•rf 

d 

T 

■cf 

WRITEPRN  [Radar  1  :=  Radar 

L  Prn  J  i 


WRITEPRN  rESM  "j  :=  ESM 
prn  J  j 
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APPENDIX  C:  SETUP  ( IFM) 


IFM  based  ESM  Receiver  case 


In  this  file  the  user  must  state  the  required  parameters  for 
both  the  IFM  based  ESM  receiver  and  the  enemy  Radar  to  be 
encountered.  The  follow  on  calculations  will  automatically 
update  their  output  every  time  a  change  is  made. 


SHIPBORNE  ENEMY  RADAR  (AIR  SURVEILLANCE) 


hr  : =  20 
Pt  :=  250000 
9 

ft  :=  3-10 

8 

310 

Xt  :=  - 

ft 

Gtm 1  :=  1585 
Gts 1  :=  .5 

PRF  :=  300 
1 

PRI  :=  - 

PRF 

-6 

PW  : =  A  ■  10 
Lt  3 
g  :=  1 
RPM  :=  6 
83dB  :=  2.5 


Enemy  radar  height  (m) 

Peak  power  (watts) 

Center  frequency  (Hz) 

Transmitter  wavelength 

Gain  Radar  Main  Lobe  (app  32  dB) 
Gain  Radar  Side  lobe  (app  -3  dB) 
Pulse  Repetition  Frequency  (Hz) 

Pulse  Repetition  Interval 

Pulse  width  (sec) 

Radar  transmission  losses 
Multipath  factor  0  <  g/'2  <  A 
Antenna  Scan  Rate 
Antenna  3dB  Beainwidth 
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60 

Ta  :=  - 

Period  of  Antenna  Scan 

RPM 

©3dB 

Time  Antenna’s  Main  Lobe  Points 

•ra  :=  Ta 

in  a  particular  direction 

360 

ESM  RECEIVER  PARAMETERS 

Vmph  :  = 

220 

ESM  platform  speed  (knots) 

Vm  :=  Vmph  ;447 

ESM  platform  speed  (m/sec) 

hesm  := 

4000 

9 

Height  of  ESM  platform  Cm) 

FL  s=  2.0  10 

9 

Minimum  Frequency  (Hz) 

FU  :=  4.0  i.O 

Upper  Frequency  (Hz) 

D  :=  FU 

-  FL 

6 

Frequency  Band  of  Coverage 

Br  s—  2000  10 

6 

RF  Bandwidth  (Hz) 

Bv  :=  10  10 

6 

Video  Bandwidth  (Hz) 

Ba  :=  2000  10 

Br 

Acceptance  Bandwidth  (Hz) 

r  :=  — 

Bv 

Ratio  of  Br  to  Bv 

gamma  : 

PG  :  = 

=  .5 

'  Br 

2 '  Bv_ 

gamma 

Integration  efficiency  coefficient 

Processing  Gain 

Fn  :=  15 

Fn 

10 

Receiver  Noise  Figure  (in  dB) 

FN  j=  10 

Receiver  Noise  Figure 
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Gesm  :=  1 

-11 

Pfa  : =  1-10 
Lesm  :=  3 
PW 

d  :  =  — 

2 

T  :=  10 
tf  :=  220 


ESM  Antenna  Gain  (Omni) 
Probability  of  False  Alarm 
Receiver  internal  losses 

Minimum  coincidence  duration 

Time  increment  (sec) 
Increment  limit 


Esm  and  Radar  input  parameters  are  written  into  output  files  called 
ESM  and  RADAR  respectively. 


i  :=  1  . . 14 


j  :  =  1  . . 20 


Radar  := 
i 


hr  ! 

pt 

ft 

xt 

Gtml 

Gtsl 

PRF 

PR  I 

PW 

R 

93dB 

Ta 

■ra 

Lt 

ESMifm  := 


!  Vmph  ! 

1  ..  I 

Vm 

hesm 

FL 

FU 

D 

Br 

Bv 

Ba 

r 

gamma 

PG 

Fn 

FN 

Gesm 

Pfa 

!  Lesm  ! 

1  '  , .  1 

1  d  j 

i  x  * 

[  —  1  t 

1  .  -  1 

I  tf  1 

UR  1TEPRN  [Radar  "j  t=  Radar 
L  Prn  J  1 

UR  1TEPRN  [ESMifm  "j  :=  ESMifm 

L  Prn  J  j 
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APPENDIX  D:  PROBABILITY  OF  DETECTION  (SUPERHET) 


Superhet  Receiver  case 

This  file  calcutates  the  probability  of-  detection  for  a 
general  ESM  System  based  on  the  theory  explained  in  Tsui’s 
book  [Ref.  33,  that  is,  considering  the  effects  of  the  ratio  Br 
to  Bv . 

Radar  :=  READPRN  [Radar  l 
L  Pm  J 

ESM  :=  READPRN  [ESM  1 

L  Prn  J 
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ESM  RECEIVER  PARAMETERS 


Vmph  : 

=  ESM 

0 

< 

3 

II 

ESM 

1 

hesm  : 

:=  ESM 

2 

FL  :  = 

ESM 

3 

FU  :  = 

ESM 

4 

D  :=  ESM 

5 

Br  :  = 

ESM 

6 

Bv  !  = 

ESM 

7 

Ba  :  = 

ESM 

8 

r  :=  ESM 

9 

gamma 

:=  ESM 

10 

PG  i  = 

ESM 

11 

Fn  :  - 

ESM 

12 

FN  :  = 

ESM 

13 

Gesm 

:=  ESM 

14 

Pfa  : 

=  ESM 

15 

Lesm 

:=  ESM 

16 

Ts  :  = 

ESM 

17 

tF  :  = 

ESM 

18 

d  :  = 

ESM 

19 

T  :  = 

ESM 

20 

tf  :  = 

ESM 

21 
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Range  Calculations 
3 

ae  :=  8493. 310  Earth  Radius  (m) 


t  :  =  0,1  . . tf 


hr 

hesm 

tl  :=  Tt 

hr  :  = 

hesm 

•  =  - 

t 

Km  1000 

J<m 

1000 

R1  : =  ae  + 

hr 

R2  :=  ae  + 

hesm 

Maximum  Range  (Radar  Horizon): 

3 


Rmax  :=  130. 34  10 


hesm  + 
Km 


hr 


Km 


(d-1 ) 


Maximum  Range  Angle: 


Brad  :=  acos 


2  2  2 
Rmax  -  R 1  -  R2 


-2R1 ' R2 


(d-2) 


6  :  = 
t 


Vm 


tl 


R1 


Range  Angle  Increment 

( d-3) 


LOS  Range  Increment: 


R  :  = 
t 


2  2 
R1  +  R2 


j2  '  R1  •  R2  '  cos  pirad  -  B 


( d-4) 


Ground  Range  increment 


R  -  hesm 
t 


( d-5 ) 
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SIGNAL  TO  NOISE  RATIO  CALCULATIONS 


-23 

k  :=  1.38-10 
To  :=  290 

Br 

B  :=  - 

eff  gamma 

f  Br  1 


2  ■  Bv 


Boltzman’s  constant 

Standard  Temperature  in 
Kelvin 

Effective  Bandwidth 

(d-6) 


N  :=  kTo-Br-FN 


Receiver  Noise  Power 

Cd-7) 


Main  Lobe  Signal  to  Noise 
Ratio : 

(at  output  of  detector) 
(d-8) 


Side  Lobe  Signal  to  Noise 
Ratio: 

(at  output  of  detector) 
(d-9) 


SNRml 

t 


2  2 

Gtml  Pt 'Gesm • g 


4 


IT 


■  N ■ Lt • Lesm 


SNRsl 

t 


SNRml  ■ 
t 


[Gts  1 


Gtml 


PROBABILITY  OF  DETECTION  CALCULATIONS 


MAIN  LOBE  SIDE  LOBE 


Kim  : =  1  +  SNRml  (d-10)  Kls  :=  1  +  SNRsl  (d-11)- 

t  t  t  t 
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SNRml  • 

t 

2 

r 

i  +  — 

4 

Main  lobe 


<d-12) 


' 

2' 

r 

1 

+  — 

SNRs 1  • 

2 

t 

2 

r 

i 

+  — 

4 

Side  lobes 


(d-13) 


- 

“ 

3  SNRml  • 
t 

•» 

2 

r 

2  +  3-  — 
4 

2 

r 

2  +  — 

4 

Main  lobe  (d-14) 


2 

r 


3  1  SNRs 1 


4 


Side  lobes  Cd-15) 


MAIN  LOBE 


SIDE  LOBES 


(d-20) 


( d-21 ) 


(d-22) 
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(d-23> 


(d-26) 


( d-27) 
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The  values  for  the  Probability  of  Detection  vs.  time  are  written 
into  output  files. 


WRITEPRN  rPdgm 

:  =  Pdm 

WRITEPRN  [*Pdgs 

r=  Pds 

Prn  . 

t 

prn 

t 
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APPENDIX  Ej  PROBABILITY  OF  DETECTION  (IFM) 


IFM  based  Receiver  case 


This  file  calcutates  the  probability  of  detection 
general  ESM  System  based  on  the  theory  explained  in  Tsui’s 
that  is,  considering  the  effects  of  the  ratio  Br  to  Bv, 


Radar  :=  READPRN  Radar 

1 


rn  J 


ESMifm  :=  READPRN  [ESMi fm 

f 


rn 


RADAR  PARAMETERS 


hr  :=  Radar 


Pt  :=  Radar 


ft  :=  Radar 


X  : =  Radar 


Gtml 

: =  Radar 

Gtsl 

:=  Radar 

PRF 

:=  Radar 

6 

PR  1 

; =  Radar 

7 

PV  : 

=  Radar 

g  :=  Radar 


iB3dB  :=  Radar 


10 


Ta  :=  Radar 


11 


ra  :=  Radar 


12 


for  a 
book. 


t  : =  Radar 


ESM  RECEIVER  PARAMETERS 


Vmph  : =  ESMi f m 


Vra  :  = 

ESMi f m 

1 

hesm 

:=  ESMifm 

2 

ii 

J 

b. 

ESMifm 

3 

FU 

ESMifm 

4 

D  :=  1 

ESMifm 

5 

Br  t  = 

ESMifm 

6 

Bv  :  = 

ESMifm 

7 

Ba  :  = 

ESMifm 

8 

r  :  = 

ESMifm 

9 

gamma 

:=  ESMifm 

PG  :  = 

ESMifm 

11 

Fn  :  = 

ESMifm 

12 

FN  :  = 

ESMifm 

13 

Gesm 

:=  ESMifm 

1 

Pfa  : 

=  ESMifm 

15 

Lesm 

:=  ESMifm 

1 

d  :  = 

ESMifm 

17 

T  :  = 

ESMifm 

16 

tf  :  = 

>  ESMifm 

19 

10 
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Range  Calculations 
3 

ae  :=  8493. 3 • 10  Earth  Radius  (m) 


t  :=  0,  1  . . tf 


hr 

hesm 

tl  :=  Tt 

hr  :  = 

hesm  :  = 

t 

Km  1000 

Km  1000 

Rl  :=  ae  + 

hr 

R2  :=  ae  + 

hesm 

Maximum  Range  (Radar  Horizon): 


3 

Rmax  :*  130. 34  10  • 

"l 

hesm  + 

hr 

:■]  Km  \ 

Km 

Maximum 


Srad  := 


Range  Angle: 


acos 


2  2 
Rmax  -  R1  - 

-2 • R1 1 R2 


(e-1 ) 


( e-2) 


Range  Angle  Increment: 


8 

t 


LOS  Range  Increment: 


tl 

t 

Vm  - 

Ri 


R 

t 


2  2 
Rl  +  R2 


2  •  Rl 


R2 ■ cos 


|8rad 


( e-3 ) 


( e-4 ) 
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SIGNAL  TO  NOISE  RATIO  CALCULATIONS 


-23 

k  :  =  1.38-10  Boltzman’s  constant 

To  :=  290  Standard  Temperature  in 

Kelvin 


Receiver  effective  Bandwidth: 


Br 

B  :=  -  <  e-5 ) 

eff  gamma 

f  Br  1 


2  ■  B  v 


Receiver  Noise  Power: 

N  :=  k  ToBr  FN  <e-6) 


Main  Lobe  Signal  to  Noise 
Ratio 

(at  output  of  detector): 


Side  Lobe  Signal  to  Noise 
Ratio 

(at  output  of  detector): 


(e-7) 


( e-8 ) 


SNRml 

t 


2  2 

G tm 1  • P t ' Gesm ' X  ' g 


4 


TT 


2 


R 

t 


•  N • Lt ■ Lesm 


SNRs  1 

t 


SNRml 

t 


[Gts  1 


|_Gtm  1 
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PROBABILITY  OF  DETECTION  CALCULATIONS 


MAIN 


Kim 

t 


K2m 

t 


K2s 

t 


K3m 

t 


LOBE 


SIDE  LOBE 


:=  1  +  SNRml  (e-9) 
t 


Kls  s=  1  +  SNRsl 
t  t 


’ 

' 

2’ 

r 

i  +  — 

1 

SNRml  ■ 

2 

1  + 

t 

2 

2 

r 

r 

i  +  — 

i  +  — 

A 

2 

' 

2 

' 

r 

l  +  — 

l 

SNRsl  • 

2 

1  + 

t 

2 

2 

r 

r 

i  +  — 

l  +  — 

_ 

. 

4 

_ 

_ 

2 


Main  lobe 


Side  lobes 


4 

3' SNRml  ■ 
t 

2 

r 

2  +  3 - 

4 

Main  lobe 

2 

r 

2  +  3  ■  — 

1  + 

2 

r 

2  +  — 

e-10) 


e-1 1 ) 


e-12) 


e-13) 


4 


MAIN  LOBE 


SIDE  LOBES 


Kim 


Kls 


K4m  :  = 


(e-15) 


K4s  :  = 


(e-16) 


t  K2m 


t  K2s 


(e-17) 


K5m  :  = 
t 


■  j  -  2  ■ ln(Pfa)  -  Kim 


(e-18) 


K5s  :  = 
t 


'•j  -  2  '  ln(Pfa)  -  Kls 


K2m 


K2s 


K3m 

1  t 


Xm  : 
t 


■•J2  ■  tr 


-•  |K4m 

“  I.  t 


*]■ 


exp 


-K4m 

max 

t 

2 

_  -500. 

(e-19) 


K2m 


K3s 

1  t 


Xs 


•j2  IT 


Ik  4  s 

L  t 


-  1  exp 


-K4s 

max 

t 

2 

L  "250 J 

(e-20) 


K2s 
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Main  lobe 


(e-21 ) 


(e-22) 


Main  lobe 


Side  lobes 


( e-23 ) 


( e-24 ) 


Bs 


Main  lobe 


<e-25) 


t 


Ys  + 
t 


2  2 


Side  lobes 


(e-25) 


Main  lobe 


Pdlm 

t 


Pd2m 

t 


Bm 

t 

Am 

t 


( e-29 ) 


Am  -  Cm 
t  t 

Am 

t 


( e-31 ) 


Bs 

t 

Pdls  :  =  - 

t  As 

t 


Pd2s  :  = 

t 


As  -  Cs 
t  t 

As 

t 


( e-30 ) 


(e-32) 


Main  lobe  probability  of  detection  (1FM  System): 

PdmIFM  :=  Pdlm  *  J.  .5  -  Pdlm  1  +  Pd2m  I  |Pd2m  -  .51  (e-33) 

t  t  L  tj  t  L  t  J 

Side  lobes  probability  of  detection  ( I FM  System): 

PdslFM  :=  Pdls  ■  s.  \ .  5  -  Pdls  1  +  Pd2s  'A|Pd2s  -  .5]  (e-34) 

t  t  L  tj  t  L  t  J 
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The  probability  of  detection  vs.  time  is  written  into  an  output 
file: 


WRITEPRN  rPdgmIFM  '  :=  PdmIFM 

prn  J  t 


WRITEPRN  rPdgslFM  ^  :=  PdsiFM 
prn  J  t 


Main  lobe 


Side  lobes 
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APPENDIX  Fj  PROBABILITY  OF  IDENTIFICATION  (SUPERHET) 


Superhet  Receiver  case 


This  file  calcutates  the  probability  of  identification  as  a 
function  of  time  Pi(t)  for  a  Superhet  ESM  System  baaed  on  the 
percent  error  deviations  of  the  input  signal  in  angle  of 
arrival,  frequency,  pulse  repetition  interval,  and  pulse  width* 
It  reads  the  Radar  and  ESM  parametrers  from  the  Radar/ESM 
parameter  file  and  creates  plots  of  Pitt)  of  both  the  emitter 
side  lobes  and  ma<n  lobe.  The  file  writes  Pi(t)  data  into  an 
output  file  for  the  final  P01<t)  calculations. 

Radar  :=  READPRN  [Radar  ] 

L  p™  J 

ESM  :=  READPRN  [ESM  ] 

L  Prn  J 


RADAR  PARAMETERS 


hr  s=  Radar 


Pt  ;=  Radar 


ft  :=  Radar 


X  :=  Radar 


Gtml  :=  Radar 


Gtsl  :=  Radar 


PRF  :=  Radar 


PR1  ;=  Radar 


PW  Radar 


g  :=  Radar 


!?3dB  :  =  Radar 


10 
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Ta 

: =  Radar 

11 

•ra 

:=  Radar 

12 

Lt 

:=  Radar 

13 

ESM  RECEIVER  PARAMETERS 


Vmph  :=  ESM 


Vm  :=  ESM 


hesm  :=  ESM 


FL  :=  ESM 


FU  :=  ESM 


D  :=  ESM 


CD 

1 

If 

ESM 

6 

Bv  :  = 

ESM 

7 

Ba  :  = 

ESil 

8 

f  ;  = 

ESM 

9 

gamma 

:=  ESM 

PG  :  = 

ESM 

11 

Fn  :  = 

ESM 

12 

FN  :  = 

ESM 

13 

Gesm 

:=  ESM 

1 

Pfa  : 

=  ESM 

15 

Lesm 

:=  ESM 

H 

Ts  :  = 

ESM 

17 

tF  :  = 

ESM 

18 

d  :  = 

ESM 

19 

10 
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T  s  =  ESM 

20 

tf  :=  ESM 

21 


Range  Calculations 

3 

ae  :=  8493. 3  10 
t  :=  0,1  . . tf 

tl  **  T- t 
t 

R1  :=  ae  +  hr 
R2  :=  ae  +  herm 

Maximum  Range  (Radar  Horizon): 


Earth  Radius  (m) 


hr 


Km 


hr 


1000 


hesm 


Rmax  : s  130 . 34 ' 10 

Maximum  Range  Angle: 


r 

- 

[hesm  + 

hr 

Km 

Km  . 

6rad  :=  acos 


2  2  2 
Rmax  -  R1  -  R2 

-2 ' R1 ’ R2 


hesm 


Km  1000 


(f-l> 


(f~2) 


Range  Angle  Increment: 
tl 


y  :  = 
t 


Vm  ' 

L  Ri 


(f-3) 


LOS  Range  Increment: 


2  2 

-1  +  R2  -  (*2 "  R1  '  R2  cos  fOrad  -  0  |1 

i  I  t  J  j 


(  f-4) 
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SIGNAL  TO  NOISE  RATIO  CALCULATIONS 


-23 

k  :=  1.23-10  Boltzman's  constant 

To  :=  290  Std.  Temperature  in  Kelvin 


Effective  Bandwidth: 

Br 

B  :=  -  (f-5) 

eff  gamma 

r  Br  1 
- 1 

2 '  BvJ 


Receiver  Noise  Power: 

N  :=  k-To  B  FN  (f-6) 

eff 


Main  Lobe  Signal  to  Noise 
Ratio 

(at  output  of  detector) 

( f  -7  > 

2  2 

Gtm 1 ' Pt  Gesm  ■  X  ■ g 

SNRml  :=  - 

t  2 

|4  n  •  R  1  N  •  Lt  ■  Lesm 


Side  Lobe  Signal  to  Noise 
Ratio 

(at  output  of  detector) 
(f-8) 


Gtsl 

SNRsl  :=  SNRml  •  - 

t  t  [Gtm  I 


Main  Lobe  Signal  Power  Side  Lobes  Signal  Power 

at  Rx. :  at  Rx . 


'  f-9) 

Sm  1  :  = 

t 


2  2 

Gtml -Pt ’Gesm  X  • g 


.  2 

4  tt’R  I  -Lt  Lesm 

tj 


(f-10) 
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PROBABILITY  OF  IDENTIFICATION  CALCULATIONS 


Frequency  error  (Main  lobe  and  Side  lobes): 


(  f  ~  12 ) 


Pulse  Repetition  interval  error  (Main  lobe  and  Side  lobes): 


dPRlm 

t 


r 

'  0.62-PRF 

:  =  I 

o 

1 

Bv  SNRml 

•1  t. 

(f-13) 
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(f-14) 


dPRIs 

t 


r 

'  0.62-PRF 

" 

j.01  -  4- 

Bv  • 

SNRsl 

. 

. 

t. 

Pulse  Width  error  (Main  lobe  and  Side  lobes): 


dPWm 

t 


r 

“ 

.7 

4 

Bv  • 

SNRml 

.01  - 

% 

t. 

. 

PW 

dPWs 

t 


.7 

4  * 

Bv  • 

SNRsl 

.01  - 

% 

t. 

PW 

- 

( f -15) 


(  f  - 16 ) 


Angle  Of  Arrival  error: 


yrms  :=  5 


Actual  rms  angular  accuracy 
(degrees ) 


Sref  :=  4 


-8 

Smin  : =  1 . 0  10 

AOA 


Angular  accuracy  of  reference 
(degrees ) 

Minimum  discernible  signal 
for  the  AOA  receiving  system 
(-80  dBw/-50  dBm) 
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(  f  - 17 ) 


(f-10) 


Identification  Matrixs 


Each  combination  is  properly  weigthed  according  to  the 
quality  of  its  components.  The  first  coefficient  is  .justified 
since  each  parameter  has  a  relative  importance,  and  therefore  it 
defines  the  overall  weight  for  each  combinations 


dAOA  =  0.3 
dF  =  0.28 
dPR I  =  0.23 
dPU  =  0.19 


And  the  second  coefficient  of  each  element  of  the  matrix 
comes  from  the  need  to  measure  relative  to  a  maximum  value  of  1 
each  of  the  combinations.  Then  if  we  just  have  obtained 
measurements  from  two  or  three  parameters,  only  one  combination 
will  be  se 1 ected . 


wl  :=  0.3  Weight  corresponding  to  AOA  error 

w2  :=  0.28  Weight  corresponding  to  freq.  error 

w3  :=  0.23  Weight  corresponding  to  PRI  error 

w4  :=  0.19  Weight  corresponding  to  PW  error 
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Main  lobe  probability  of  identification  matrix: 


(  f-19 


Pirn 

t 


max 


1 ' 

fo.  25  ■  fdAOAm  +  dFm  +  dPRlm  + 

dPWm  ]] 

L  L  t 

t 

t 

tjj 

-.fwT 

+  w2  +  w3 '  fo.333' 

dAOAm  + 

dFm 

+  dPRlm  11 

L 

t 

t 

tjj 

■jwl 

+  w2  +  w4  ‘  fo .  333  • 

dAOAm  + 

dFm 

+  dPWm  'll 

I 

L. 

t 

t 

t  < 

•jwl  +  w3  +  w4 '  [b.333'  dAOAm  +  dPRlm  +  dPWm 

L  L  t  t  tj. 


\|w2  +  w3  +  w4-  fb .  3333 '  fdFm  +  dPRlm  +  dPWm  1] 

L  L  t  t  tjj 


■■jwl  +  w2  ■  fb.  5  '  fdAOAm 

+  dFm  11 

l  L  » 

tjj 

■>jwl  +  w3  •  JO.  5  •  JdAOAm 

+  dPRlm  Jj 

■•jwl  +  w4  ‘  fo.  5  fdAOAm 

+  dPWm  11 

L  L  t 

tJJ 

■■jw2  +  w3  ‘  p.5-  JdFm  + 

dPRlm  11 
tjj 

■Jw2  +  w4  ■  fo.5-  fdFm  + 

dPWm  11 

L  L  t 

tjj 

•Jw3  +  w4  ■  fo .  5  ■  fdPRIm 

+  dPWm  11 

L  L  t 

tjj 

••jw2  1  dFm 

t 

■jw3  -  dPRlm 

t 


r — 

•|w4  ■  dPWm 


t 
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Side  lobes  probability  of  identification  matrix 


( f -  20 ) 


Pis 

t 


=  max 


1-  [b.  25  '  fdAOAs  +  dFs  +  dPRIs  + 

L  L  t  t  t 

d™\]] 

>  w  1  +  wC  +  w3-  fo.  333  1  fdAOAs  +  dFs 

+  dPRIs  I! 

1 - 

C* 

r* 

tJJ 

•■Jwl  +  w2  +  w4  ■  fo.333-  fdAOAs  +  dFs 

+  dPWs  1] 

L  L  t  t 

tJJ 

-•jwl  +  w3  +  w4  •  fo.333  -  fdAOAs  +  dPRIs 

+  dPWs  1] 

L  L  t 

t  tJJ 

■jw2  +  w3  +  w4-  jo.  3333-  jdFs  +  dPRIs  +  dPWs 


’jwl  +  w2  •  fo.  5  ■  fdAOAs 

+  dFs  11 

L  L  t 

tJJ 

'jwl  +  w3  ‘  fo .  5  •  fdAOAs 

+  dPRIs  11 

L  L  t 

tJJ 

>wl  +  w4  1  fo.  5  ■  fdAOAs 

+  dPWs  11 

1 - 

C+ 

tJJ 

■jw2  +  w3  '  fo.  5  ’  fdFs  + 

dPRIs  11 

L  L  t 

tJJ 

•jw2  +  w4  ■  fo.  5  ■  fdFs  + 

dPWs  11 

L  L  t 

tJJ 

>jw3  +  w4  ■  fo.  5  ’  fdPRI  s 

+  dPWs  11 

L  L  t 

tJJ 

■fw2dFs 

t 

■jw3  ■  dPR  I  s 

t 

-i: 


w4 • dPWs 


t 
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The  Probability 
an  output  file; 


of  identification  vs 


time 


is  written  into 


WRITEPRN 

Pi  dm 

:=  Pirn 

WRITEPRN  [P ids 

:=  Pis 

prn 

t 

.  ’  Prn  J 

t 
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APPENDIX  G;  PROBABILITY  OF  IDENTIFICATION  ( IFM) 


IFM  based  Receiver  case 


This  file  calcutates  the  probability  of  identification  as  a 
function  of  time  Pi(t)  for  an  IFM  -  ESM  System  based  on  the 
percent  error  deviations  of  the  input  signal  in  angle  of 
arrival,  frequency,  pulse  repetition  interval,  and  pulse  width. 

It  reads  the  Radar  and  ESM  parametrers  from  the  Radar/ESM 
parameter  file  and  creates  plots  of  Pi(t)  of  both  the  emitter 
sidelobes  and  mainlobe.  The  file  writes  Pi(t)  data  into  an 
output  file  for  POI(t)  calculations. 

Radar  :=  READPRN  [Radar  1 

L  Prn  J 

ESMifm  :=  READPRN jESMi fm  ] 

L  prn  J 


RADAR  PARAMETERS 


hr  : 
Pt  : 
ft  : 
X  :  = 
Gtm  1 
Gts  1 
PRF 
PR  I 
PU  : 


=  Radar 

0 


=  Radar 

1 

=  Radar 

2 


Radar 

3 

:=  Radar 


4 

:=  Radar 

5 

:=  Radar 

6 

:=  Radar 

7 

=  Radar 

8 


g  :=  Radar 

S 

93dB  :=  Radar 

10 


112 


Ta 

:=  Radar 

11 

■ra 

:  =  Radar 

12 

Lt 

:=  Radar 

13 

ESM  RECEIVER  PARAMETERS 


Vmph  :=  ESMifm 


Vm  :=  ESMifm 

1 

hesm  :=  ESMifr 
FL  :=  ESMifm 

3 

FU  :=  ESMifm 


D  :=  ESMifm 


Br  s=  ESMifr 
Bv  ;  =  ESMifm 

i 

Ba  :=  ESMifm 
f  J=  ESMifm 


6 


8 


gamma  :=  ESMifm 


PG  :=  ESMifm 

I 

Fn  :=  ESMifm 


10 


11 


12 


FN  :=  ESMifm 

13 

Gesm  :=  ESMifm 

14 

Pfa  :=  ESMifm 

15 

Lesm  :=  ESMifm 

16 

d  :=  ESMifm 

17 

T  :=  ESMifm 

18 

tf  ESMifm 

19 
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Range  Calculations 


3 

ae  :=  8493. 3'10  Earth  Radius  (ml 


t  : =  0,1  . . tf 


hr 

hesm 

tl  :=  T  t 

hr  :  = 

hesm  := 

t 

Km  1000 

Km  1000 

R1  :=  ae  + 

hr 

R2  :=  ae  + 

hesm 

Maximum  Range  (Radar  Horizon): 


3 

Rmax  :  =  130. 34  10  • 

" 

hesm  + 

hr 

Km 

Km 

Maximum 


6rad  := 


Range  Angle: 

2  2 
Rmax  -  R1 


acos 


-2R1 ■ R2 


( g-1  > 


(g-2) 


Range 


0 


Angle  Incremen 

f 


I  Vm  - - 


t : 


(g-3) 


t  L  R1  J 

LOS  Range  Increment: 


R  :  = 

2  2 

R 1  +  R2  -  J 

2  R1  R2  cos  fftrad  - 

8  n 

( g~4 ) 

t 

L 

tjj 
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SIGNAL  TO  NOISE  RATIO  CALCULATIONS 


-23 

k  : —  1.3810 
To  :=  290 


Effective  Bandwidth: 

Br 

B  := - 

eff  gamma 

T  Br  1 


2  •  Bv 


Receiver  Noise  Power: 


N  :=  k-To-B  • FN 
eff 


Boltzman’s  constant 

Std.  Temperature  in  Kelvin 


(g-5) 


(g-6) 


Main  Lobe  Signal  to  Noise 
Ratio 

(at  output  of  detector) 
(g-7) 


Side  Lobe  Signal  to  Noise 
Ratio 

(at  output  of  detector) 
(g-8) 


2  2 

Gtml  Pt'GesmX  g 

SNRml  :=  - 

t  2 

j4  ■  tt  ■  R  j  ■  N  ■  Lt  •  Lesm 


SNRsI  :=  SNRml 
t  t 


Gtsl] 

Gtm  I  J 


Signal  Power  at  ESM  from 
Radar  mainlobe  : 


Signal  Power  at  ESM  from 
Radar  Side  lobes  : 


( g-9 ) 


Sm! 


2  2 

Gtml  Pt 'Gesm  A  1 g 


4  tt-R 

t. 


•  Lt • Lesm 


(g-10) 


Ss  1 

t 


Sm  I 

t 


Gts  1 
Gtml .. 
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PROBABILITY  OF  IDENTIFICATION  CALCULATIONS 


6 

Frequency  Resolution  =  2  MHs.  ; < 1 1  bit  word  length)  R  :=  210 
Frequency  error  (Main  lobe  and  Side  lobes)  : 


dFm  :=  2 
t 


■ 

• 

■ 

r  1 

A  , 

"1 3 

4 

max 

SNRml  ■ 

2  3 

TT  -PW  B 

j 

t 

ef  f 

.001  - 

- 

R 

- 

ft 

. 

(g-ll> 


dFs  :=  s. 
t 


.001  - 


max 


-17 


2  3 

SNRsl  -tt  -PW  B 

t  ef  f 

R 


ft 


( g- 12 ) 


Pulse  Repetition  interval  error  (Main  lobe  and  Side  lobes): 


dPRlm 

t 


:=  2 

■ 

.01  -  4- 

"  0. 62  PRF 

j 

1 

Bv •  SNRml 

’1 

(g-13) 
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dPRls  5  .01  -  4 

t 


0.62  -PRF 


Bv  SNRsl 
•1  t 


(g-14) 


Pulse  Width  error  (Main  lobe  and  Side  lobes) 

r  Lf _ :7  m 


dPWm  :=  5  .001  - 
t 


Bv  ■ 

SNRml 

tj 

(g-15) 


.7  1" 

Bv-  SNRsl 


dPWs  ;=  I  .001  - 
t 


(g-16) 


Angle  Of  Arrival  error 
STms  :  =  5 

Sref  :=  4 


Smin  :=  l.o  10 
AOA 


Actual  rms  angular  accuracy 
(degrees ) 

Angular  reference  accuracy 
(degrees ) 

Minimum  discernible  signal 
for  the  AOA  receiving  system 
(-80  dBw  or  -50  dBm) 
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I 

dAOAm  : = 

4.0 

•  2  fSm  1 

-  Smin  i 

<  g-17> 

t 

9rms_ 

L  t 

AOA  J 

dAQAs 

t 


4.0 

-  -5  [Ssl 

©rms  J  L  t 


Smin 

AOA 


(g-18) 


Identification  Matrix: 


Each  combination  is  properly  weigthed  according  to  the 
quality  of  its  components.  The  first  coefficient  is  justified 
since  each  parameter  has  a  relative  importance,  and  therefore  it 
defines  the  overall  weight  for  each  combination: 


dAOA  =0.3 
dF  =  0.26 
dPR 1  =  0.23 
dPU  =  0, 19 


And  the  second  coefficient  of  each  element  of  the  matrix 
comes  from  the  need  to  measure  relative  to  a  maximum  value  of  1 
each  of  the  combinations.  Then  if  we  just  have  obtained 
measurements  from  two  or  three  parameters,  only  one  combination 
will  be  se 1 ected . 


wl 

ii 

O 

0) 

Weight 

correspond i ng 

to 

AOA  error 

w2 

:=  0.28 

Weight 

correspond i ng 

to 

Freq.  error 

w3 

:=  0.23 

Weight 

corresponding 

to 

PR1  error 

w4 

:  =  0.19 

Weight 

corresponding 

to 

PW  error 
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Main  lobe  Probability  of  identification  matrix: 


( g~ 19 ) 


Pim  :=  max 
t 


1-  (b .  25  ■  IdAOAm  +  dFm  +  dPRlm  +  dPUm  11 

L  L  t  t  f  tJJ 

•■jwl  +  w2  +  w3-  jo.  333-  JdAOAm  +  dFm  +  dPRlm  JJ 

Jwl  +  w2  +  w4-  jo.  333-  JdAOAm  +  dFm  +  dPWm  JJ 

■jwl  +  w3  +  w4-  jo .  333  •  jdAOAm  +  dPRlm  +  dPWm  JJ 
Jw2  +  w3  +  w4-  jo.  333-  jdFm  +  dPRlm  +  dPWm  JJ 
wl  +  w2 '  To.  5  ■  fdAOAm  +  dFm  ]] 

L  L  t  tJJ 

•jwl  +■  w3-  jO.  5-  jdAOAm  +  dPRlm  JJ 
■Jwl  +  w4  •  jO.  5-  JdAOAm  +  dPWm  JJ 
•jw2  +  w3  ■  jO.  5  •  jdFm  +  dPRlm  JJ 

•^w2  +  w4 '  [0.  5  ■  [dFm  +  dPWm  "]~j 

L  L  t  tJJ 

Jw3  +  w4  ■  jo .  5  ■  jdPRlm  +  dPWm  J'j 
•jw2  ‘dFm 


dPR  I  m 


dPWm 
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Side  lobes  Probability  of  identification  matrix: 


( g-20 ) 


f  1  ’  [O.  25  •  fdAOAs  +  dFs  +  dPRls  +  dPUs  T 

L  L  t  t  t  tj. 


Pis  :=  max 

t 


•jwl 

+  w2  +  w3  • 

0.333-  fdAOAs 

+  dFs 

+  dPRls  1] 

L  L  t 

t 

tJJ 

'|wl 

+  w2  +  w4 • 

0.333-  fdAOAs 

+  dFs 

+  dPWs  I! 

L  L  t 

t 

tJJ 

I  w  i 

+  w3  +  w4  1  fc 

1.333-  fdAOAs 

+  dPRls 

+  dPWs  11 

L 

L  t 

t 

tJJ 

-■Jw2 

+  w3  +  w4  ■ 

0.333-  fdFs  + 

dPRls 

+  dPWs  11 

L  L  t 

t 

tJJ 

•jwl 

+  w2 

'|wl  + 

w3 " 

'|wl  + 

w4  • 

•|w2 

+  w3 

■•jw2 

+  w4 

•  w3  + 

w4  ■ 

JdAOAs  +  dPRls 
0.5  ■  [dAOAs  +  dPWs 

L  t  t 


J] 


■fw2 ' 

■|w3  ■ 
'|w4  • 


dFs 


dPR  I  s 


dPWs 
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Z5jP£°mi!Ity  °f  identifi<:ation  vs-  tin.  is  written  into 


WRITEPRN  fPidmlFM 

:=  Pim 

WRITEPRN  [Pids I FM  1 

L  prn  J 

t 

L.  prn  J 

APPENDIX  H:  PROBABILITY  OF  COINCIDENCE  (SUPERHET) 


Superhet  Receiver  case 


This  file  calcutates  the  probability  of  coincidence  for 
a  Superhet  ESM  -  Frequency  hopper  (pulse  to  pulse)  Radar  type  of 
encounter  based  on  the  theory  explained  in  Wiley's  book,  that  is 
considering  an  approach  using  window  functions. 

Radar  :=  READPRN  [Radar  1 

L  pm  J 

ESM  :=  READPRN  [ESM  1 

L  p™  J 


RADAR  PARAMETERS 


hr  :=  Radar 

0 

Pt  :=  Radar 

1 

ft  :=  Radar 

2 

A  : =  Radar 

3 


Gtm  1 

: =  Radar 

Gts  1 

:=  Radar 

PRF 

:=  Radar 

6 

PRI 

:=  Radar 

7 

PW  : 

=  Radar 

8 


g  :=  Radar 

9 

03dB  :=  Radar 


10 

Ta  :=  Radar 

11 

•ra  :=  Radar 

12 

Lt  :=  Radar 


13 
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ESM  RECEIVER  PARAMETERS 


Vmph  :=  ESM 


Vra  :=  ESM 


hesm  :=  ESM 


FL  :=  ESM 


FU  :=  ESM 


D  :=  ESM 


Br  :=  ESM 


Bv  :=  ESM 


Ba  :=  ESM 


T  :=  ESM 


gamma  : *  ESM 

10 

PG  :=  ESM 

11 

Fn  :=  ESM 

12 

FN  :=  ESM 

13 

Gesm  :=  ESM 

14 

Pfa  :=  ESM 

15 

Lesm  :=  ESM 

16 

Ts  :=  ESM 

17 

Tf  :=  ESM 

18 

d  :=  ESM 

19 

T  :=  ESM 

20 

tf  :=  ESM 

21 
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Range  Calculations 


3 

ae  :=  8493. 3'10  Earth  Radius  (m) 


t  :  =  0, 1  .  .  tf 

hr  hesm 


tl  :=  Tt 

hr  :  * 

hesm  := 

t 

Km  1000 

Km  1000 

R1  :=  ae  + 

hr 

R2  :=  ae  + 

hesm 

Maximum  Range  (Radar  Horizon): 


3 

Rmax  :=  130.3410  • 

hesm  + 

hr 

\  Km  • 

Km 

Maximum  Range  Angle: 


6rad  :=  acos 


2  2  2 
Rmax  -  R 1  -  R2 


•2  R1 - R2 


( h-2 ) 


Range  Angle  Increment: 


0 

t 


LOS  Range  Increment: 


2  2 

R  :  = 

R1  +  R2  - 

2  R1  ■  R2  •  cos  ferad  -  £ 

1 

t 

!  ! 

,  L  t. 

J 

(h-3) 


( h-4 ) 
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PROBABILITY  OF  COINCIDENCE  CALCULATIONS 


n  : 

=  5 

Number 

•rl 

:  =  ra  -  d 

( h-6 ) 

■r2 

:  =  rf  -  d 

<  h-B ) 

t3 

:=  PW  -  d 

o 

H 

1 

JC 

of  frequencies  hopped  over 

t4  :=  2 ■ PW  -  d 

■rli  :=  (Ta  -  ra)  -  d 

Thop  :=  n’PRI 


Mean  time  between  coincidences: 

Ta  Ts 'PRI -Thop 

Tom  :=  - 

t!-t2-t3  +  t! ■ t2 ■ fA  +  t1-t3"t4  +  r2r3  r4 


Tos 


Ta  Ts  PR!  Thop 

■rl  1  •  r2  •  r3  +  rll  r2  t4  +  •ril,  r3-T4  +  t2,t3-t4 


Instantaneous  Coincidence  Probability: 


■rl  ■  r2  •  t3  ■  t4 

Pom  :=  - 

TaTs-PRI -Thop 


•rl  1  '  r2  ■  r3  •  r4 

Pos  :=  - 

Ta  Ts -PRI -Thop 


Km  :=  1  -  Pom  <h-16)  Ks  : =  1  -  Pos 


(h-5) 

(h-7) 

( h-9  > 

( h-1 1 ) 


( h-12  > 

( h-13) 


(h-14) 

(h-15) 


( h-17) 
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Probability  of  Coincidence  vs.  Time: 


[-“111 

t  l  l 


Pcml  :=  1  -  Km-exp  max  - 

t  Tom 

L  L-600JJ. 


(h-18) 


Pcs  1 

t 


|Ks -exp  max 


-tl  IT 
t 

Tos 

-6OOJJ. 


(h-19) 


Probability  of  coincidence  vs.  time  Is  written  into  an  output 
file: 


WR1TEPRN  jPcmhop  1  ;=  pcm  1 

L  Prn  J  t 


WRITEPRN 


Pcshop  1  :=  Pcs  1 
prn  J  t 
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APPENDIX 


PROBABILIY 


OF  COINCIDENCE  (IFM) 


IFM  Receiver  case 


This  file  calcutates  the 
based  ESM-Radar  type  of 
in  Wiley’s  book,  and  in 
considering  an  approach 


probability  of  col 
encounter  based  on 
Self's  article  [Ref 
using  window  funoti 


ncidence  fo 
the  theory 
s.  2  and  12] 
one. 


r  an 
expla 
that 


IFM 

ined 

is, 


Radar  :=  READPRN  [Radar  1 

L  Prn  J 

ESMifm  s=  READPRN  [ESMifm 

L  pm 


RADAR 

PARAMETERS 

hr  :  = 

Radar 

0 

Pt  :  = 

Radar 

1 

*-* 

II 

Radar 

2 

X  :  = 

Radar 

3 

Gtm  1 

:=  Radar 

4 

Gts  1 

:=  Radar 

5 

PRF  : 

=  Radar 

6 

PR  1  : 

=  Radar 

7 

PW  :  = 

Radar 

8 

g  !  = 

Radar 

9 

83dB 

:=  Radar 

10 

Ta  :  = 

Radar 

11 

•ra  *  = 

Radar 

12 

Lt  :  = 

Radar 

13 
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ESM  RECEIVER  PARAMETERS 


Vmph  :=  ESMifm 


Vm  :  = 

ESMifm 

1 

hesm 

:=  ESMifm 

2 

FL  :  = 

ESMifm 

3 

FU  :  = 

ESMifm 

4 

D  :  = 

ESMifm 

5 

Br  :  = 

ESMifm 

6 

CD 

< 

II 

ESMifm 

7 

Ba  :  = 

ESMifm 

8 

r  :  = 

ESMifm 

9 

gamma 

:=  ESMifm 

PG  :  = 

ESMifm 

11 

Fn  :  = 

ESMifm 

12 

FN  :  = 

ESMifm 

13 

Gesm 

:=  ESMifm 

1 

Pfa  j 

=  ESMifm 

15 

Lesm 

:=  ESMifm 

1 

d  ;  = 

ESMifm 

17 

T  :  = 

ESMifm 

18 

tf  :  = 

ESMifm 

19 

10 


Range  Calculations 


3 

ae  :=  8493. 310  Earth  Radius  (a) 


t  :  =  0,  1  . . tf 


hr 

hesm 

tl  :=  Tt 

hr  :  = 

hesm  :a 

t 

Km  1000 

Km  1000 

R1  : =  ae  + 

hr 

R2  :=  ae  + 

hesm 

Maximum  Range  (Radar  Horizon): 


0 

o 

Rmax  :=  130.3410  • 

hesm  + 

hr 

Km 

Km 

Maximum  Range  Angle: 


Brad  :=  acos 


2  2  2 
Rmax  -  R1  -  R2 

-2R1 ' R2 


(i-2) 


Range  Angle  Increment: 
tl 


W  :  = 
t 


Vm 


R1 


( i-3) 


LOS  Range  Increment; 


R  :  = 
t 


2  2 
R1  +  R2 


J2  ’  R1  ‘  R2  cos  jSrad  -  6  j"j 


( i-4) 
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PROBABILITY  OF  COINCIDENCE  CALCULATIONS 


n  J=  5  Number  of  frequencies  to  be  hopped  over  (i-5) 


■rl 

:  =  ra  -  d 

(i-6) 

•t3  :=  2  •  PW  -  d 

(  i 

-7) 

•r2 

:=  PW  -  d 

(i-8) 

Thop  :=  nPRI 

(i 

-9) 

rll  :=  (Ta  -  1 0 ■ ra )  -  d  (i-10) 


Mean  time  between  coincidences: 


Tom 


Ta  ■  Thop ■ PR  I 


•rl  r2  +  rl  r3  +  r2  •  r3 


Main  lobe  (i-11) 


Tos 


Ta-Thop-PRI 

•rl  1  ■  t2  +  rll  r3  +  r2  •  t3 


Side  lobe 


( i  -12) 


Instantaneous  Coincidence  Probability: 


•rl  ■  r2  •  r3 

Pom  :=  -  Main  lobe  (i-13) 

Ta  Thop  PR  1 


'rll  r2.  t3 

Pos  := -  Side  lobe  (i-14) 

Ta -Thop -PR  I 


Km  : =  1  -  Pom 


(i-15)  Ks  :=  1  -  Pos  (i-16) 
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Probability  of  Coincidence  vs.  Time: 


Pcm 1 IFM 

t 


Mai n 1 obs 


Pcs  1 IFM 

t 


Side  1 obe 


(i-17) 


(  i  -18) 


Probability  of  coincidence  vs.  time  is  written  into  an  output 
file: 


WR1TEPRN  iPcml IFM  n  :=  Pcm 1 IFM 

L  prn  J  t 


WRITEPRN  (Pcsl IFM  1  :=  Pcs  1 IFM 
prn  J  t 
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APPENDIX  J;  PROBABILITY  OF  INTERCEPT 


Thi9  file  reads  all  other  probability  files  and  calculates  the 
final  Probability  of  Intercept  as  function  of  time,  PQltt),  fflT 
various  conditions. 


GET  THE  ITERATION  TIME  AND  STEP  SIZEl 

ESM  :  ==  READPRN  |ESM  ] 

L  prn  J 

ESMifm  :  =  READPRN  [ESMi fm  ] 

L  prn  J 


T  :=  ESM  tf  :=  ESM 

20  21 

t  :=  0, 1  .  .  tf 

tl  :=  t-T 
t 


1.-  READ  ALL  THE  OTHER  FILES 

a.)  Superhet  files: 

Pdgm  :=  READPRN  jPdgm  "j 

L  Prn  j 

Pdgs  :  = 

READPRN  jPdgs  ] 

L  Prn  J 

Pcmhop  :=  READPRN jPcmhop  ~j 

L  Prn  J 

Pcshop 

:=  READPRN iPcshop 

L  P 

Pidm  :=  READPRN  [Pi dm  "j 

L  Prn  j 

Pids  := 

READPRN  j'P ids  ] 

L  Prn  J 
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b. )  I FM  files: 
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b. )  Probability  of  Intercept  ( Si de 1 otaes )  1 


PO I  1 f ms  :  =  Pdgg 1 FM  -PcslIFM  -PidsIFM 
t  t  t  t 


(  j-5) 


c, )  Overall  Probability  of  Intercept  (IFH  Receiver): 


POIifm  :=  PO  1 1  f  mm  +  POIifms  -  fPOUfmni  -POIifms  ]  (j-6> 
t  t  t  L  t  t. 
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